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Abstract 

We analyze the thermodynamic properties of interfaces in the three-dimensional Falicov 
Kimball model, which can be viewed as a primitive quantum lattice model of crystalline 
matter. In the strong coupling limit, the ionic subsystem of this model is governed by 
the Hamiltonian of an effective classical spin model whose leading part is the Ising Hamil- 
tonian. We prove that the 100 interface in this model, at half-filling, is rigid, as in the 
three-dimensional Ising model. However, despite the above similarities with the Ising model, 
the thermodynamic properties of its 111 interface are very different. We prove that even 
though this interface is expected to be unstable for the Ising model, it is stable for the Falicov 
Kimball model at sufficiently low temperatures. This rigidity results from a phenomenon of 
"ground state selection" and is a consequence of the Fermi statistics of the electrons in the 
model. 
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1 Introduction 

Domain boundaries can have important effects on the transport properties of condensed 
matter materials. In some cases, transport is believed to occur mainly or exclusively along 



domain boundaries ||35|| . This is related to the drastic effect the presence of a domain wall, or 
interface, can have on the low-lying excitations of the system. Therefore, it is important to 
study non-periodic equilibrium states, in particular interface states, in statistical mechanics, 
and to understand their stability, fluctuations, and low-lying excitations. 

Many phenomena in condensed matter physics, which are of current interest, are intrinsi- 
cally quantum mechanical in origin. Quantum effects also play a crucial role in the properties 
of interfaces. There are, however, very few rigorous results on interface states in quantum 
statistical mechanics. In [§, |3| a general perturbation theory was developed, which, under 
certain assumptions, shows that a small quantum perturbation does not destroy an existing 
interface Gibbs state of a classical discrete spin model (the so-called Dobrushin states). 

In this paper we consider the opposite situation, namely, one in which quantum fluc- 
tuations stabilize the interface against thermal fluctuations, while the classical limit does 
not have a stable interface. We demonstrate the occurrence of such an "order by disorder" 
effect J|3L |2lfl , in a simple quantum lattice model — the three-dimensional Falicov Kimball 
(FK) model [See Section |2| for a description of the model]. We are motivated by recent work 
P5| which demonstrates that, at zero temperature, in the two-dimensional ferromagnetic 
XXZ Heisenberg model with Ising-like anisotropy the quantum fluctuations stabilize the 11- 
interface (i.e., an interface in the diagonal direction). In the present work we prove such an 
effect at finite temperature for two interfaces in the FK model. 

The FK model was chosen because its statistical mechanical properties have been studied 
extensively (see, e.g., [19|, |29|, |22], |20], [L8|, [33], [27]) and, recently, convenient perturbation 



expansions have been developed for it f33|, [XT], |23|, p2| . We expect that the XXZ model in 



three dimensions also has a stable 111 interface at sufficiently low temperatures. Its analysis 
is, however, more involved due to the presence of gapless excitations in the interface |26|, [3l| . 

Our main result is that, in three dimensions, the FK model has a stable 111 interface 
at sufficiently low temperatures. This should be compared with the three-dimensional Ising 
model since, in the strong coupling limit, the FK model can be considered as a perturbation 



of the Ising model (see e.g. [29, 33, Q). Dobrushin showed that the Ising model has 100- 
interface states, but its Ill-interfaces are expected to be unstable at any finite temperature. 
It has been proved recently that in the zero-temperature limit of the three-dimensional Ising 
model the 111 interface fluctuates p5| . This is related to the degeneracy of the ground states 
with a Ill-interface which grows exponentially with the volume (the rate of exponential 
growth may depend on the boundary conditions however! See ^TJ). In the FK model this 
degeneracy is lifted. In this sense, this is an example of the phenomenon of "ground state 
selection" [^l|] by quantum fluctuations. We refer to Section |27| for a detailed discussion of 



the 111 interface configurations. 



In [24] Lieb and Kennedy showed how to study the FK model in terms of an Ising- type 



model for the Ising configurations that is obtained by taking the trace over the electron states 
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for any given ion configuration. The Hamiltonian for the ions can be explicitly computed to 
any order in perturbation theory together with a bound on the sum of the higher order terms 
[see Appendix A and ||32|| ] . For the study of the 100-interface one needs the explicit form 
of this Hamiltonian up to second order. For the Ill-interface fourth order perturbation 
theory is needed. In principle, our method could be used to study interfaces with more 
general orientations, but higher order terms in the perturbation series will be needed; e.g., 
the 112 interface is infinitely degenerate at fourth order, but we expect it to be stabilized 
at sufficiently low temperatures by the sixth order terms. Therefore, one should expect that 
the Falicov-Kimball model has an infinite number of interface phase transitions. 

We follow Dobrushin |TJ| (see also |5|, |6|]) in proving the existence of an interface by con- 
sidering an effective two-dimensional model for the interface. Although, this two-dimensional 
model turns out to be quite complicated, and involves many-body interactions of arbitrar- 
ily long range, we can analyze it with a Peierls-type argument. It is probably possible to 
develop a general Pirogov-Sinai theory P§L EO, ESI H to treat this situation, along the lines 



of [|12j and f3~7 ], but we found it more convenient to make efficient use of simpler methods. 
The result is a more transparent and relatively short proof. 

Our main results are stated at the end of the next section. Our main technical tool is 
the convergence of certain cluster expansions proved in Appendix B. Appendix A contains 
the proof of a bound on the remainder term in the expansion of the effective Hamiltonians 
for the ions. In Section ^ we discuss the rigidity of a 100-interface, which is much simpler 
than the case of the Ill-interface treated in Section [|. 

2 The FK model and effective Hamiltonians 



The Falicov Kimball (FK) model |TB| is a lattice model of spin-polarized electrons and 
classical particles (ions) [16| . The electrons and ions interact via a purely on-site interaction. 
The electrons can hop between nearest neighbour sites of the lattice, but the ions are static. 
There is a hard-core repulsion between the ions, which prevents more than one ion from 
occupying a single lattice site. The presence or absence of an ion at a lattice site x is 
described by a classical variable, W(x), which is equal to unity if there is an ion at the site 
x and is zero otherwise. Let the number operator of an electron at the site x be given by 
n = c x c x , where c x and c x are the creation and annihilation operators of the electron at the 
site x. Let /i, and /i e denote the chemical potentials of the ions and electrons respectively. 
Let Z 3 denote an infinite cubic lattice, with unit lattice spacing, such that the coordinates 
of the sites are given by half-integers. For notational simplicity in our description of the 
interface, it is more convenient to consider this lattice instead of Z 3 . The Hamiltonian of 
the model defined on a finite lattice A C Z 3 is given by 

n A (t,U)=H 0A (U)+tV A , (2.1) 



where 1~Coa{U) is a Hamiltonian which is given entirely in terms of an on-site interaction as 
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follows: 

zGA xeA x£k 

:= E $ o, (2.2) 

The operator V\ causes electrons to hop between nearest neighbour sites of the lattice: 

Va = -t ( c l c v + c l c x), 

<xy>cA 

E V x , (2.3) 

X=<xy>cA 

Here < xy > denotes a pair of nearest neighbour sites in the lattice. The hopping amplitude 
of the electrons is denoted by t G 1R. The first rigorous study of the above model was done by 



Kennedy and Lieb ||24|| . They considered the classical particles to be nuclei and the on-site 
interaction to be the Coulomb attraction between nuclei and electrons. In accordance with 
this interpretation, they chose the coupling constant U to be negative. They proved that 
the ground states of this model display perfect crystalline ordering for the choice 

IH = = U, (2.4) 



the ions being arranged in a checkerboard configuration. The choice Q2.4D corresponds to 
a neutral model: the average number of electrons in the lattice is equal to the the average 
number of ions, both being equal to half the number of lattice sites (half-filling). 

The model described by the Hamiltonian T~C\(t, U) for U < 0, is mathematically equiva- 
lent to the model with U > (see and there is a simple relation between the properties 
of the attractive model, U < 0, and the repulsive one, U > 0. In this paper we work with 
U>0. 

We study the FK model in the strong-coupling limit i.e., for U » \t\, and hence consider 
the hopping term tV\ to be a perturbation to the Hamiltonian TCoa(U)- It is convenient to 
renormalize the hopping amplitude t to unity. This amounts to the following rescaling: 

U/t — ► U 
(3t — (3, (2.5) 

where (3 = l/ksT (&b being the Boltzmann constant and T the absolute temperature). In 
our expansions U^ 1 plays the role of a small parameter. 

For a fixed value of the coupling constant U, the zero temperature phase diagram of the 
unperturbed Hamiltonian 7{$\(U) in the plane of chemical potentials can be easily obtained 
|l8| . To obtain the ground states of Hqa(U), for any given set of values the chemical 



potentials /ij and /z e , we only need to find the single-site configuration which minimizes 
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$ x (|2.2|) . When both the chemical potentials are negative, the ground state corresponds 
to all sites of the lattice being empty. In the rest of the (fa- fa,) plane, it is found that for 
values of the chemical potentials such that fa, < 2U and/or fa < 2U, there is no doubly 
occupied site at zero temperature. For these values of the chemical potentials, the ground 
state corresponds to an all-ion configuration if fa > fa. with fa > 0, and to an all-electron 
configuration if fa, > fa with fa, > 0. For < fa = fa, < 2U all singly occupied configurations 
are equally likely and hence the ground state is infinitely degenerate. The origin and the 
point fa = fa, = 2U also correspond to infinitely many ground states. At the origin, each site 
is either empty or singly occupied, whereas at the point (2U, 2U) each site is either singly- 
or doubly occupied. At zero temperature, for fa > 2U and fa, > 2U, every site is doubly 
occupied by an electron-ion pair. 

As in [^, [J5J, I52I , we will rely on the fact that for U > c\t\, where c is a positive constant, 
the ionic subsystem of the FK model defined on a finite cubic lattice A can be described by 
an effective classical Hamiltonian. We will study only the neutral model at half- filling,, i.e., 
fJ'i — fJ'e — U. Then, it follows form the circuit representation of |3j| and the bounds proved 



in Appendix A that the equilibrium states of the ions are described by an effective classical 
Hamiltonian of the form 

nf(U) = ±- £ s' x s' y + K^,U), (2.6) 

<xy>eA 



where 



s' x := 2W(x) - 1. (2.7) 



The variable s' x can be interpreted as an on-site spin variable since s' x = 1 if there is an ion 
at x and s' x — — 1 otherwise. Note that TIa(U) depends on f3 through the remainder term 
TZ\(f3, U). This is unavoidable if one wants an exact correspondence between the correlation 
functions of the ions in the Falicov-Kimball model and the same correlation functions of the 
effective classical spin model. In particular, if one writes TZ\(/3, U) as a sum of products of the 
s' x variables, one sees that it still contains a nearest neighbour contribution with a coefficient 
that tends to zero exponentially fast as (3 — > oo. As all our results are for f3 and U sufficiently 
large, this temperature dependence will be of no consequence however. Similarly, as we will 
do later on [( gjg ) and (|2~35D l, one can extract from TZ\(P,U) the leading contributions 



for the next-nearest neighbour and plaquette interactions, which are independent of (5, and 
estimate the temperature dependent corrections with the bounds proved in Appendix A. 

All terms of higher orders in the perturbation parameter U , i.e., all terms of order 
U~ n , with n > 3, are contained in the remainder TZ^{f3, U), which, just as the leading term, 
depends on (3 in an inessential way. Hence we simply write it as TZ\(U). This remainder is 
expressible in terms of local classical interactions {Rb(U)}: 

Ka(U)= ]T R b (U) (2.8) 

|fl|>2 
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Here B denotes a connected set of lattice sites, i.e., if x, y £ B then there exists a sequence 
of sites Xq — x, Xx, . . . , x n — y such that Xi G B and \xi — Xi + \\ — 1 for all i — 0, . . . , n — 1. 
The number of sites in B is denoted by \B\. We refer to such a set B as a bond. A bond can 
be represented by a graph, the vertices being the sites and the lines of the graph representing 
nearest neighbour bonds between pairs adjacent sites. Let B denote the set of all bonds in 
the lattice A. For each bond B appearing in the above sum (|2.8| ), the interaction Rb{U) can 
be expressed as a product of two or more on-site spin variables s' x with x G B. 

Note that, for us, a bond is, by definition, a connected set. This is convenient in com- 
binatorial arguments and a natural choice in view of the way perturbation theory produces 
long-range and multi-body interactions as a composition of nearest neighbour hoppings. It 
does not exclude the presence of terms of the form s x s' y , with \x — y\ > 2, in the effective 
Hamiltonian. Such terms are included in connected bonds B containing x and y. 

A configuration on Z 3 , denoted by lu, is therefore given by a set of assignments {s' x } xe z 3 
of s' x G { — 1,1} to each x G Z 3 . For any finite subset Y C Z 3 , let ujy denote the restriction 
of the configuration uj to the subset Y . A boundary condition (b.c.) will be specified by a 
configuration uj, meaning that for any finite volume A G Z 3 , the system is considered with 
a fixed value for the spins on each x G Z 3 \ A, determined by uj. 

From Lemma [T9| of Appendix A (in particular (|A.75|) ) it follows that for (3 and U large 
enough, there exist positive constants, c\ and C2 (with C\/U < 1) such that: 

• for \B\ > 3 

W.G0I < *QT B) (2-9) 

where 

g(B) :=n(B)-l, (2.10) 

with n(B) being defined as the minimum length (in units of lattice spacing) of a closed 
path which passes through all sites of B. 

• while, for \B\ — 2, 

\Rb(U)\ < ~C2(^)\ (2.11) 
The latter bound [( |2.11| )1 results from the following fact: We have that g{B) = 1 for \B\ = 2; 



however, the term in the effective Hamiltonian of order Z7 -1 , has been extracted and is given 
by the first term on the RHS of (^^). The contribution from electron hoppings between 
nearest neighbour sites, to the remainder term, are therefore of order U~ 3 , as even powers 
of U~ l do not occur in the expansion. 

The above bounds imply that there exists a constant r > such that 

J2 \R B (U)\e r9{B) < oo (2.12) 

B30 
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and hence the interaction {Rb(U)} decays exponentially. Using the bound (|2.11|) , the effec- 
tive Hamiltonian (p. 61) can be written as follows 



nf(u) = n^(u) + nl 3 (u) (2.13) 

with 

n$(u) = -L £ 44+ £ R B (j3,u) 

<xy>cA snA^0 
|fl|=2 

= m e 44 ( 2 - 14 ) 

<xy>eA 

where 

J(U):=J(P,U)^±j + h.o.>0, (2.15) 

the symbol "h.o." denoting terms of higher orders in U^ 1 which are bounded by all" 3 , for 
some constant a > 0. Throughout the rest of the paper, whenever the symbol h.o. appears 
in a sum involving powers of U" 1 , it will denote the presence of a correction term which is 
bounded by a positive constant times the next odd power of U~ x . The correction term is an 
infinite series that can be computed order by order in perturbation theory, and that has a 
dependence on (3. This temperature dependence is inessential and we will routinely omit it 
from the notations. We will only use that, for /3U sufficiently large, a bound of the form ( |2.9| ) 
holds. We refer to Appendix A for the proofs of the bounds. A more general class of models 
as well as a detailed discussion of the temperature dependence of the effective Hamiltonian 
is contained in . 



We define 

Kf(U):= £ R B {U) (2.16) 



|S|>2 



(2) 

The leading part T^oav^O °f the effective Hamiltonian is identical to the Hamiltonian of an 
antiferromagnetic Ising model with nearest neighbour interactions of strength J(U) in the 
presence of a magnetic field of strength h. 

The fact that the Hamiltonian 7i A s (U) ( [2.13D is invariant under "spin-flip", is a conse- 
quence of half-filling. The leading part, H^(U) [(gjj)], of the effect ive Hamiltonian has 



two ground states in each of which the +'s and — 's occupy alternate sites of the lattice, 
i.e., the antiferromagnetic Neel states. Thus, to order U~ l , the effective Hamiltonian gov- 
erning the ionic subsystem has a two-fold degenerate ground state, while for the original 
unperturbed Hamiltonian Hoa(U) (given by (|2.2j )) the ground state energy is independent 
of the ion configuration. Hence, as far as the ionic subsystem is concerned, the effect of 
the quantum perturbation tV is to select two ground state configurations from the infinitely 
many ground states of H.oa(U). This phenomenon in which the quantum perturbation lifts 
the infinite degeneracy of the classical ground states is known as "ground state selection" 
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zjj . Moreover, for low enough temperatures, the characteristic long range order of the 
Neel states persists under the action of the remainder IZ^ 3 |)TT |. The main purpose of this 



paper is to prove that a similar phenomenon of ground state selection occurs for the 111 
interface of the FK model and that the selected interfaces are rigid, i.e., the interfaces persist 
in the thermodynamic limit at finite but non-zero temperature. See Section [5]. 
It will be more convenient for us to perform the transformation 

*x— ►s* = (-l) 2N < fora11 x=(x 1 ,x 2 ,x 3 )eA, (2.17) 

where \x\ = x% + x<i + x%. The above transformation yields an equivalent Hamiltonian with 
the leading part given by the ferromagnetic Ising model: 

Ht™{U) = -J{U) s * s y> (2-18) 

<xy>eA 

where J(U) is given by (|2.15|) . Let us denote s x = 1 by the symbol "+" and s x = —1 by the 
symbol "— ". The two ground states of Ti { ^ ro {U) are denoted by s + and s_ and correspond 
to each site in the lattice being occupied by a + and a — respectively. We refer to these two 
ground state configurations and their low-temperature analogues as homogeneous phases. 
For a finite volume Ac2 3 , the boundary conditions defined 

either by s x = +1 for all x E Z 3 \ A, 

orbys x = -1 for all x e Z 3 \ A. (2.19) 

are referred to as homogeneous b.c. 

It is an interesting question what boundary configuration of the effective spin system 
correspond to bona fide boundary conditions for the original Falicov-Kimball model without 
introducing special boundary interaction terms. It is not hard to see from the derivation 
of the effective Hamiltonian in p2| , that the homogeneous b.c. discussed above, as well as 
the boundary conditions employed later in this paper to construct interface Gibbs states, 
can be achieved by imposing Dirichlet boundary conditions for the electrons on a volume 
with includes the first boundary layer. On the effective Hamiltonian, this has the effect 
of truncating the interactions across the boundary to nearest neighbour interactions only, 
which is inconsequential. This simple correspondence between boundary conditions, however, 
is only possible because we are considering spin-less fermions, at half-filling and with nearest 
neighbour hopping only. 

In more general situations natural boundary conditions in the original system will lead to 
modified boundary interactions in the effective spin system. The general strategy adopted 
in 



p2| to address this problem is to trace over all possible configurations of static fermions 
outside the volume A under consideration. Due to the Pauli principle, each such configura- 
tion of spinless fermions simply defines an excluded volume. This means that the effective 
Hamiltonian describes a weighted average of Gibbs states obtained with different Dirichlet- 
type boundary conditions for the fermions. As the fermions can wander across the boundary 



Rigidity of FK Interfaces 



9 




Figure 1: The intersection of a finite cubic lattice Ac2 3 centered at the origin, 
with a plane passing through the origin and orthogonal to the vector n = (1, 1, 1). 



of A, the effective Hamiltonians thus obtained will have interaction terms of arbitrary range 
across the boundary. As is clear from our analysis here, this kind of "averaged" bound- 
ary conditions can equally well be used to demonstrated to existence of Gibbs states with 
interfaces. 

As mentioned before, the purpose of this paper is to study the thermodynamic properties 
of the 100 and 111 interfaces of the three-dimensional FK model. In finite volume, the 
presence of an interface can be enforced by a suitable choice of mixed boundary conditions 
in the standard way ]T^] (see eqs. fl2.20| ) and ( |2.21| ) below). 

To construct the 100 interface Gibbs state we consider A C Z 3 to be a parallelepiped 
centered at the origin. The 100 interface is orthogonal to the vector n = (0,0,1). The 
boundary condition which leads to a 100 interface (which we shall refer to as "b.c. 1") is 
given as the configuration {s x } on the sites x — (xi, x 2 , x 3 , ) G A c := Z 3 \ A denned by 

a =( +1 if *3>V2 (2 20) 

Sx l-l ifx 3 <-l/2 [Z - Z()) 

For an analysis of the 111 interface consider the intersection of a a plane passing through 
the origin and orthogonal to the vector n = (1, 1, 1) with a cube A C Z 3 . This intersection 
yields a plane bounded by a hexagon (as shown in Figure p]) which divides A into two equal 



volumes. To obtain a 111 interface we consider the spin variables s x [ (|2.17|) 1 to have opposite 
values on the two sides of this dividing plane. 

More precisely, the boundary condition that leads to a 111 interface (which we shall refer 
to as "b.c. 2") is given as the configuration {s x } on the sites x = (xi,X2, £3, ) G A c = Z 3 \ A 
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defined by 

+1 if Xl + x 2 + x 3 > 1/2 



1 if xi +x 2 + x 3 < -1/2 i 2 - 21 ) 

Each of these boundary conditions divide the volume A into two subvolumes (the config- 
uration in the latter being given by the two homogeneous phases s + and s_) and enforce 
the occurrence of an interface (domain wall) pinned to the boundary of the volume. The 
residual free energy per unit area of this interface is its surface tension and is denoted by the 
symbol r mixcdb c -. it is defined by taking the difference of the free energy of the system in 
the volume A under mixed b.c. and the free energy corresponding to the homogeneous b.c. 

1 " mixed 

r mixedb.c. = u _ l 0g (=A_) ( 2 .22) 

Ayzs P\1a\ k C a ' 

where 5™ lxed and S^ om are the partition functions of the system in the finite volume A w.r.t 
mixed b.c. and homogeneous b.c. respectively (see Section |3|). The symbol |/a| denotes the 
area of the portion of the ground state interface which is contained in the volume A. Such 
a definition is justified since the volume contributions proportional to the free energies 



of the coexisting phases, as well as the boundary effects, cancel and only the contributions 
to the free energy of the interface are left. 

The limit limAyz 3 is taken in a definite order : e.g. if the ground state interface lies 
in a plane P n , passing through the origin, which is orthogonal to a non-zero vector n, then 
the dimensions of the lattice perpendicular to the plane P n are taken to infinity before the 
dimensions parallel to it are taken to infinity. The symbol |/a| denotes the area of the portion 
of the interface which is contained in the volume A. 

It will be convenient to consider the relative Hamiltonian defined with respect to the 
homogeneous phases s + and s_. For any configuration s = {s x } x( zz3 the relative Hamiltonian 
is given by 



H A (s) := -J(U) £ (s x s y -l)+ £ $fl(a) 

<xy>eA snA^0 

S|>3 

:= H^(s)+Rf(s) (2.23) 



where 



$ B (» := R B (s) - R B (s+) = R B (s) - Rb{s~), (2.24) 

and J(U) is given by ( |2.15|) . Here and henceforth the explicit [/-dependences of the relative 
Hamiltonian and its components have been suppressed for notational simplicity. 

Each configuration in a finite volume A, with respect to fixed boundary conditions (ho- 
mogeneous or mixed) can be geometrically described by specifying the Ising contours which 
are defined as follows p4| : We define a face to be a unit square which bisects a nearest neigh- 
bour bond of the lattice perpendicularly. To each nearest neighbour bond we can associate a 
face. A face / belongs to A if at least one site of the corresponding nearest neighbour bond 
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is in A. Given a configuration uj\ on a finite lattice A, with b.c. u, let S w (ua) be the set of 
faces associated with nearest neighbour bonds between opposite spins. Decompose ^(oja) 
into maximally connected pairwise disjoint components. Each such component is referred 
to as an Ising contour (or simply contour, if confusion is not likely) and is denoted by the 
symbol 7. For homogeneous b.c. the contours are closed surfaces lying entirely within the 
volume A. However, for mixed b.c. there is necessarily one (and only one) contour which 
is pinned to the boundary of the volume A. This is the only infinite maximally connected 
component of the set iS w (c<ja) and is referred to as the interface. 

There is a one-to-one correspondence between spin configurations on the lattice and non- 
intersecting families of contours T = {7}: uo A = ^(r). We shall refer to such families as 
compatible families of contours. In the sequel we shall use the symbol 7 to denote both the 
contour and its support. The number of faces in a contour 7 is denoted by I7I and satisfies 
the bound I7I > 6, since the lattice is three-dimensional. The energy of a contour 7 is given 
by 

^(7):=^A ) (r / U{ 7 })-< ) (r'). (2.25) 

for any set, T', of non-intersecting contours, not containing 7, such that T'IJ{7} is again a 
family of non-intersecting contours. 

It follows from ( [2.23| ) that the relative energy of a configuration u>a(T) is given by 

H A (T) = H$(T) + Rf(T) (2.26) 

with 



where 



< ) (r) = E^(7), (2-27) 

7gr 



£(7) = 2J(U) = (^ + h.o.)| 7 | =: Ji(U)\l\ (2-28) 



-2U 
and 

i?| 3 (r) = £ $ B (r). (2.29) 



BnlV0 

|S|>3 

The condition B R V 7^ denotes that the above sum runs over all bonds B such that 
B n 7 7^ for some 7 G T, i.e., the bond B intersects at least one face of a contour 7 G T. 

The quantity 22(7) is the self-energy of the contour 7 w.r.t. H^. Its definition (|2.28|) implies 

(2) 

that, for sufficiently large U, the Hamiltonian ifg A satisfies the Peierls condition: 

E(ri) > , (2.30) 



with a Peierls constant 



J = c U-\ 



(2-31) 
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where cq is a positive constant. Moreover, from the definition fl2.24| ) and the bound (|2] 
follows that for \B\ > 3 

. (Ci\9(B) 



it 



1$ 



B\ 



(2.32) 



with C2 = 2c2. The effect of the term Rj^{U) is to modify the self-energy of the contours and 
also to introduce interactions between the contours. Hence the spin model is reformulated 
as a model of interacting contours. 

This contour Hamiltonian can be used to prove the rigidity of the 100 interface. However 



we show in Section |2~3"1 that, under the boundary condition b.c.2 (|2.21|) , the leading part f/. 



OA 

of the Hamiltonian yields infinitely many ground state interfaces in the thermodynamic limit. 
These interfaces are characterized by the fact that they all have minimal area. We refer to 
such interfaces as minimal area interfaces. Hence, to prove the rigidity of the 111 interface 
we consider a more detailed decomposition of the relative Hamiltonian in which all the terms 
up to order U~ 3 of the perturbation series are computed explicitly and retained in its leading 
part H${: 



rr(4) , d>5 
-"OA + -"A 



(2.33) 



where 



H, 



(4) 
OA 



-fi- 



ll 



-hi,,,) Yl < 

<xy>cA 



U£7 16U 3 



S X Sy 



K S X Sy 



\x-y\=2 



16U 3 



+ h.o.) (s x s y s z s t -l) (2.34) 

x,y,z,tcP(A) 



where -P(A) is the set of plaquettes, each plaquette consisting of four lattice sites forming a 
unit square. 

The remainder R^ 5 is obtained from the series ( p.29|) defining R^ 3 by subtracting all 
terms which depend on U~ n with n < 3. It is given by 



R 



>5 



E * 

|S|>3 



Si 



(2.35) 



with the potentials <3># satisfying the bound 



1$ 



B 



< c 2 10 



rn(B) 



[2.36) 



where 



m(B) := max(5, g{B)). 



[2.37) 



In the expression for the relative Hamiltonian Hq2 ( p.23|) the terms up to order U 1 



arc 



computed explicitly by second order perturbation theory, while the terms of order U in 
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Hq2 ( |2.34|) are obtained by fourth order perturbation theory. We shall refer to (|2.23|) as 
the second order decomposition of the relative Hamiltonian, and ( |2.33| ) as its fourth order 
decomposition. 

In Section [| we prove that from the infinite set of minimal area interfaces, and up to 
translations, a unique interface configuration is selected (i.e., attributed minimal energy) 
by Hq . This interface and its translations are referred to as ground state interfaces of 
the three-dimensional FK model under the boundary condition b.c.2 ( |2.21|) . Hence the 111 
interface exhibits ground state selection. Further, we prove that the selected interface is 
rigid in the sense that it persists under the action of the remainder i?^ 5 (£7) at sufficiently 
low temperatures, in the thermodynamic limit. 

Our main results are that, for sufficiently large U, and at sufficiently low temperatures, 
the Gibbs states obtained in the thermodynamic limit A /* Z 3 with the boundary conditions 
b.c. 1 and b.c. 2, describe rigid interfaces in the 100- and 111 directions respectively. For 
a precise statement of our main results we introduce the following notations: Let <>[b.c.i] 
and <> [b.c.2] denote the expectation values in the (infinite-volume) Gibbs states with the 
mixed boundary condition b.c.l (|2.20|) and b.c.2 ( |2.21| ) respectively. Further, we recall that 
for a site x = (xi,x 2 ,x 3 ) in the lattice, s x denotes the on-site spin variable defined through 
and ( |2.17|) . Using these notations and definitions, we state our main results through 



the following theorems: 

Theorem 1 There exist positive constants Uq and D such that for all U > Uq, and (3/U > 
D , the following bounds are satisfied: 

< s x >[b. c .i]> 1 " 2C e- c 'P /u for x 3 > 1/2, (2.38) 

and 

< s x > [b . c .i]< -1 + 2Coe- c,f3/u for x 3 < -l/2„ (2.39) 
where Co and d are positive constants given in terms of Uq and D . 



Theorem 2 There exist positive constants Uq, D and D'q such that for all U > Uq, (3/U > 
Dq and fl/U 3 > D'o, the following bounds are satisfied: 

< s x > [b .c.2]> 1 - 2{C Q e- c " (3 ' u3 + C ie - c "V u } for x x + x 2 + x 3 > 1/2, (2.40) 

and 

< s x > [b . c .2]< -1 + 2{C Q e- c " (3 ' u3 + C 1 e-^' u } for Xl + x 2 + x 3 < -1/2, (2.41) 
where Co, c" , C\ and c[ are positive constants given in terms ofUo, Do an d D'o- 
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To prove these results we follow the method introduced by Dobrushin [15] and consider 



effective two-dimensional models of the 100- and 111 interfaces, obtained by projecting the 
interfaces on the planes defined by x 3 = and x\ + x 2 + x 3 = respectively [see Section [I] and 
|2.2| for details] The rigidity of the interfaces follows from an analysis of the low-temperature 
properties of these effective two-dimensional models. 

2.1 The geometry of 100 interfaces. 

The geometry of the 100 interfaces of the FK model is the same as the geometry of 
the 100 interfaces described by Dobrushin for the three dimensional Ising model. Hence we 



refer to [15] for the definitions of geometrical objects which describe the interfaces and their 
significances as configurations of a two dimensional (contour) model: the ceilings (which 
project on to ground states of the two dimensional model), the walls (which project on to 
the contours), and the standard walls (which project on to the external contours). 

2.2 The geometry of 111 interfaces. 

The geometry of the 111 interfaces of the FK model is much more involved. Let X denote 
the family of interfaces under the b.c.2 ( |2.21| ) and I denote its typical element. Such an 
interface is pinned at the boundary of A on the curve defined by 

<9AH {(x,y,z) G Z 3 | x + y + z = 0} (2.42) 

In this section we describe the underlying geometrical structure necessary for the definition 
and study of the effective two-dimensional model for such an interface. 

For each integer n we define planes P„ C Z 3 C R 3 orthogonal to the vector n = (1, 1, 1), 
by P n = {(x, y, z) 6 Z 3 | x + y + z = n}. Let P denote the plane in 1R 3 which contains Pq. 
Let V denote the orthogonal projection onto P and let Pa denote the portion of the plane 
P which is contained in the volume A. 

An effective two-dimensional model for the 111 interface is obtained by an orthogonal 
projection of the interface / onto P. Its complete description requires the following ingredi- 
ents: 

1. A set of vertices V = Vo U V\ U V2, where V nmo d3 = V(P n ), 11 G Z, is a triangular 
lattice in P with lattice constant v2. The set V also forms a triangular lattice, but 



with lattice constant y 2/3. 

2. A set of edges £ = £01 U £12 U £20, where S{j is the set of nearest neighbour edges in 
the lattice Vj U Vj. All edges have Euclidean length y^2/3. 

3. A set of triangles T consisting of all elementary triangles in V. 



4. 



A set of rhombi 1Z = IZo U TZ± U 7^2, where TZi is the set of all rhombi formed by two 
triangles in T that share an edge e G £{oi2}\{i}, * = 0,1,2. 
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Figure 2: A pair of rhombi in good position. The edge shared by the two rhombi 
is referred to as a good edge. 



The set of vertices V is the projection of the vertices in 2 3 , and V((x,y, z)) E V% if and 
only if i = (x + y + z) mod 3. 

The set of edges £ is the projection of the set of nearest neighbour bonds in H? . They 
are the edges of a triangular lattice with lattice constant v/2/3. For each pair of distinct 
i, j G {0,1,2}, the set of vertices V« U Vj forms a triangular lattice TCij (also with lattice 
constant yj2/3), and with edges Together the three cover V twice. 

As before, to each nearest neighbour bond in the lattice Z 3 we associate a unit square 
(face) which bisects it perpendicularly. Recall that the vertices of such a face have integer 
coordinates. The rhombi in TZ are the projections of these faces. Hence, an interface / in 
the lattice A projects onto a covering of the plane Pa with rhombi in TZ. We refer to such 
a rhombus covering as a rhombus configuration or, for brevity, as an ^-configuration (to be 
distinguished from a configuration of Ising contours). 

• A rhombus is said to be an overlapping one if it contains the projection of more than 
one face of the interface. Otherwise it is said to be non-overlapping. Each triangle t in 
an ^-configuration C\ necessarily belongs to the projection of an odd number of faces 
of the interface. To each such triangle t we associate a number o(t) which we refer to 
as its overlap number, and define as follows: 

o(t) := {the number of faces of the interface whose projection contains t} — 1 (2.43) 

A triangle with a non-zero overlap number is referred to as an overlapping triangle. It 
is evident that each overlapping triangle has an even overlap number. 

• Two non-overlapping rhombi which share an edge are said to form a good pair if the 
angle enclosed by them is 27r/3, i.e., if up to translations and rotations the pair is as 
shown in Figure 0. The edge shared by such a pair is referred to as a good edge. We 
consider the good pairs as open complexes. This means that a good pair is composed 
of an open edge together with the two adjacent open rhombi. Two good pairs are 
connected if their intersection is an open rhombus. 
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• We define the type of a rhombus r, r(r), to be j if r G 7^, with j G {0, 1, 2}. 

Let a ta/mg of the plane P be defined as a complete covering of P with non-overlapping 
rhombi in TZ. For each pair of distinct i,j G {0,1,2} the set of edges S \ Sij drawn in 
the plane P yield a tiling of P with the rhombi in TZ{ox2}\{ij} ■ Minimal area interfaces and 
ground state interfaces have simple geometric descriptions in terms of tilings (see Section 

The mixed boundary condition b.c.2 translates into a boundary condition for the R- 
configuration in Pa- It is given by a tiling of P \ Pa with rhombi of a single type, say TZq- 
We call this the standard b.c. for the rhombus model. 



2.3 The minimal area interfaces 

An interface / is of minimal area if and only if its projection V(I) is a tiling of P. 
Equivalently each such tiling is in one-to-one correspondence with a dimer covering of the 
hexagonal lattice dual to the triangular lattice V, i.e., of the lattice with set of sites given 
by the centers of the triangles in T. It is less obvious, although also a well-known fact in 



enumerative combinatorics [£8], that all tilings of P with rhombi in TZ correspond to 
a unique minimal area interface. This is shown in the following proposition. The proof is 
constructive, i.e., it provides an algorithm for obtaining the interface from the tiling and vice 
versa. 

Proposition 3 The tilings of the plane P with rhombi in TZ under standard b.c. are in 
one-to-one correspondence with the minimal area interfaces in the volume A under the mixed 



boundary condition b.c.2 ^2.2f\ ). 



Proof: As noted above it is obvious that each minimal area interface under the boundary 
condition b.c.2 ( f2.21| ) projects onto a tiling of Pa with rhombi in TZ. Hence, we only need 
to show that to each tiling there corresponds exactly one interface that has that tiling as 
its projection. The interface will automatically be minimal. This amounts to associating a 
unique face of Z 3 to each rhombus in the tiling such that the resulting set of faces form a 
connected set which is pinned at the boundary of A along the curve defined by Q2.42j ). The 



projection of the set of faces constituting an interface / yields a set of rhombi in TZ that 
covers Pa. For each face / we can number its vertices ax, a 2 , 03, in such a way that there 
is a unique integer n(f) for which 

ax G P n (/)_x, a 2 G P n (f), «3 G P„(/) + i, a 4 G P n (f) 

It is easy to see that if V(f) = r G TZi, i — 0, 1, 2, then i = n(f) mod 3, and that V(f) and 
n(f) uniquely determine /. 

For any tiling with standard boundary conditions we will construct a unique height func- 
tion h : V — > Z with the property that for each rhombus r in the tiling the heights of its 
vertices, when ordered appropriately, and such that {vi,Vi + x} are edges of r, are given by 



h(vx) = n — 1, h(v2) = n, h^v^) — n + 1, h(v±) = n 
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Figure 3: The labeling of the six unit vectors emanating from a site of the triangular 
lattice. 



for some integer n satisfying % — n mod 3 iff r G TZi. It follows that h satisfies \h(v) — h(w)\ = 
1 for each edge {v,w} of a rhombus in the tiling. A minimal area interface in A under the 
b.c.2 ( |2.21|) is an interface whose projection on the plane Pa is a tiling. It consists of all 
faces {/} in A for which: 

1. V(f) is a rhombus in the tiling, and 

2. the vertices v { of V(f) satisfy {h( Vi ) | 1 < i < 4} = {n(f) - 1, n(f), n(f) + 1}. 

It remains to construct the height function h and to verify that it is the unique function 
with the stated properties. Let us denote by e\, . . . , e 6 the vectors of minimal length (= \J\), 
emanating from a single point of the triangular lattice (spanned by the edges in S) such that 
the tips of these vectors are the vertices of a hexagon. See Figure ||[ For any two vertices 
v, w G V for which {v, w} is an edge in a tiling of P, w — v is one of e«. Let us denote by 
#(e*j, Cj) the angle between e*j to e}, which is a multiple of n/3. We claim that for each tiling 
there is a unique height function h satisfying 

h( \ h( \ — / ^ 5 w — v) is an even multiple of ir/3 
\ — 1 if 9(e\, w — v) is an odd multiple of n/3 

for each edge {v, w} of a rhombus in the tiling. The height function h must then be obtained 
by summing up the above differences along edges, starting from a convenient reference value 
on the boundary of Pa- Consistency of this definition follows from the elementary observa- 
tions that 

1. any two paths connecting the same pair of vertices and consisting of edges in the tiling, 
together enclose a bounded subset of the plane tiled with rhombi, and 

2. for any rhombus the differences h(w) — h(v) along the four edges of a rhombi sum up 
to zero because the rhombus has two angles of 7r/3 and two of 27r/3. 

Then the uniqueness of h is also obvious. 

Note that the height function h has been defined such that the height of any vertex in 
the tiling is equal to the sum of the coordinates of the point in the interface of which it is 
the projection. | 
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Proposition 4 The number of minimal interfaces Na grows exponentially with the area of 
the interface and satisfies the bounds: 

2 1 -¥<N A < 2 2 ^ 
where |/a| is the area of the minimal area interfaces in the volume A. 

Proof: Consider the interface defined by the TZo tiling, and pick one of the three hexagonal 
sublattices of the vertices. Note that, independently of each other, each hexagon can be tiled 
in two ways with three rhombi. This proves the lower bound. 

The upper bound is obtained by considering the hexagonal lattice dual to the triangu- 
lar lattice. Since every tiling is in one-to-one correspondence with a dimer covering of this 
hexagonal lattice, it is also in one-to-one correspondence with a path covering of the hexag- 
onal lattice, restricted by the condition that at each site the incidence of the path is two. 
The upper bound is just the straightforward bound on the number of paths given by 

2 l where / is the length of the path 

As the path is covering, the length of the path equals the number of bonds in the hexagonal 
lattice which equals twice the number of rhombi needed to tile the region, i.e., twice the 
area. H 

Note that the exact rate of exponential growth of the degeneracy depends on the shape 
of the finite volumes [[11]. 



From Proposition |] it follows that in the limit A /* Z z there are infinitely many minimal 
area interfaces. 

3 The relevant partition functions 

We shall define the partition functions of the FK model in terms of the second- and 



fourth order decompositions of the relative Hamiltonian. These are defined through ( 2.23 ) 
- ( |2.24| ) and ( |2.33p - ( |2.35| ) respectively. Such definitions are justified because the quantity 
relevant for the study of an interface is not a solitary partition function but rather the 
quotient of two partition functions, namely, a partition function corresponding to mixed b.c. 



and one corresponding to homogeneous b.c (see ( 2.22 )). The use of relative Hamiltonians 
in the definition of the partition functions corresponds to the simultaneous subtraction of 
the energy of a homogeneous phase from the Hamiltonians appearing in the numerator and 
denominator of the quotient and hence keeps the quotient unchanged. 

Let and be the partition functions in A C 2 3 w.r.t the homogeneous boundary 
conditions defined by (|2.19|) . The spin-flip symmetry of the effective Hamiltonian 7i A (£7) 
( |2.13| ), for the choice h = 0, implies that 



^+ ^— ^hom (n ^ 
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Let S A 00 and H A n be the corresponding partition functions under the mixed boundary con- 
ditions b.c.l ( |2.20|) and b.c.2 ( [2.21| ) respectively. 

A mixed b.c. (b.c.l ( |2.20| ) or b.c.2 ( |2.21| )) leads to the appearance of an interface / 
which divides the volume A into two subvolumes A" and A b j which are, respectively, the 
regions above and below /. The configuration in these subvolumes are defined by finite sets 
of compatible contours T a : = {7", ...,7^} and T b = {7J, ...,7^} respectively. We define 
r := T a U T b and denote a configuration on the lattice by (JUT). Let X and X denote the 
families of interfaces resulting from b.c.l and b.c.2 respectively. Let / and / denote their 
typical elements. 

In terms of the second order decomposition ( [2.26] ) of the relative Hamiltonian H\, the 
partition function for homogeneous b.c. is given by 



5A m = E n e ~ /3Jl(C/)l71 ' II e-^ s(r) , (3.2) 
r= {7l ,..., 7 „}cA i=i w» 

where J\{U) ~ ^ + h.o. The partition function relevant for the mixed boundary condition 
b.c.l ( p.20| ) is also obtained by using eqns.( |2T26| ) - ( |2.29| ). It is given by 

2a° = E E e~^ Iur \ (3.3) 
iex r : =r a ur b 

where Si(I U T) is the energy of the configuration (/ U T). It is the value that the relative 
Hamiltonian H\ takes on the configuration (JUT): 

^(/ur)= E Mu)\Y\ + E Ji(u)\i b \ + Ji(u)\i\+ E s*(/ur) (3.4) 

7 a gr a 7 i 'gr 6 B:\B\>3 

sn(/ur)#0 



The first three terms on the r.h.s. of ( |3.4| ) are the energies of the contours in the configuration 
(/ U T) as defined through (|2.25|) . The third term is the energy of interaction among these 



contours and arises from the long-range tail potential R A ( |2.29p of the relative Hamiltonian 
H A . 

Let (/ U 0) denote a configuration which consists only of the interface / and no other 
contours. The energy of such a configuration can be interpreted as the "bare" energy of the 
interface, i.e., the energy of the interface in the absence of any other contour. Let us denote 
this energy by E^™^). It is given by 

E^ IC (I):=J 1 (U)\I\ + E <M/U0). (3.5) 

B:\B\>3 
Bnl^fb 



From ( |2.28p and ( |2.32j ) it follows that for sufficiently large U, the energy i?J 3are (7) satisfies 



the bound 

\E^ e (I)\ < const. ^ | /| (3.6) 
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It is convenient to isolate the "bare" energy of the interface from the remaining terms in 



the expression (|3.4j) for 8\{I U Y). From (|3.4j) and ( |3.5|) it follows that 

£i(I U T) = Y E(j) + E b ^ c (I) + E 1 (IU r), (3.7) 

7er 



with 



^i(Jur) = Y <M/ur)- £ $ B (/u0) 



fl:|S|>3 B:\B\>3 

sn(iur)/0 sn/^J 



]T {$ B (/ur)- x (Bn/^)$ B (/u0)} 



B:\B\>3 

sn(iur)^0 



:= £ $' B (/Ur), (3.8) 

B:|B|>3 

fln(/ur)^0 

where %(•) denotes the characteristic function. Note that 

$^(/U0) = O if 5nr = 0. (3.9) 

Hence, the functions Q' B for which B intersects only the interface, do not contribute to the 
energy Ei(lUT). The contributions of such bonds is included in the "bare" energy E hare (I) 
of the interface. 

A non-zero $' B (/U0) arises only from those bonds B which intersect at least one contour 
in T. This observation allows us to write E\(I U T) as follows. 

E 1 (IUT)= Y WUT) (3.10) 

B:\B\>3 

Further, it follows from ( |2.32p that <&'b(I U T) satisfies the bound 



|$' B (/ur)| < |$ B (/ur)| + |$ B (/u0)| 

<Cl\9(B) 

-u 



< (3.11) 



for \B\ > 3. For bonds B which do not intersect the interface /, 

$s(/ur) = $ B (/ur) (3.12) 

The partition function for the boundary condition b.c.l, defined by ([Of) is hence given 



by 

: ioo = ( .;/.' " /; x ST Y ' i '"' un 



"A 

i& r"={7f,...,7#}cA| r 6 ={ 7l b ,...,7 g 6 }cAj 



f[ e -vmW x [[ e -«l (3.13) 

3=1 fc=l 
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with Ei are (I) and Ei(lUT) being denned through (|375|) and ( p.lQ ) respectively. The partition 
function for the boundary condition b.c.2 (|2.21| ) can be similarly written as 

Sl n = E E e~^\ (3.14) 

where £ 2 (I U T) is the energy of the configuration (JUT). In order to determine whether 
the 111 interface is rigid, it is necessary to use the fourth order decomposition [eqns. (|2.33| ) 
- ( [2.35 )] of the relative Hamiltonian H\, for computing the contribution of the interface / 
to the energy £2(1 U T). However, it is sufficient to consider the second order decomposition 
[eqns. ( [2.261 ) - (|2.29|) 1 for evaluating the corresponding contribution of the contours in T. We 
first introduce some notations which are useful in evaluating £ 2 {I U T). 

• Let Pj C -P(A) denote the set of plaquettes in Z 3 which are intersected by the faces 
of /. Let p denote a typical element in this set. 

• Let Bj C Z 3 denote the set consisting of pairs of next nearest neighbour sites, 

{{x,z} E Z 3 \ \x-z\ = 2}, (3.15) 
such that the line joining each pair is intersected by a face in /. 
From ( |2.33|) - (|2.35| ), it follows that the energy £ 2 {l U T) is given by 

£ 2 (IUT) = j a (r/)|j| + (_L + h.o.) £ ^ } (j) + (_L_ + h. .) £ h p (I) 

{x,z}£Bj p={x,y,z,t}£Pj 

+ £ l B (/ur)+ Yl Ji(u)\r\+ E -Mt0l7T + E ^(/ur), 

(3.16) 



B:\B\>3 -yagpa 7 b gr i> B:\B\>2 

sn;>0 " snr/0 

sn/=0 



where 



h p := 5(s x s y s z s t - 1) + 3(s x s z + s y s t - 2) (3.17) 
for each p = {x, y, z, t} G P/, 

h{ x , z } := s x s z - 1 (3.18) 

for each {x, z} e and 

Let (J U 0) denote a configuration which has the interface I as its only contour. We denote 
the corresponding "bare" energy of the interface by E 2 are (I). It is defined as follows. 

£ 2 bare (/) := J 2 ([/)|J| + (-L+h.o.) £ h {x , z} {I) 

{x,z}(i)eBj 

+ ( T^ + h ' a) £ W+ E $b(/U0). (3.20) 

p^X^tyePf B:\BJ>3 
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The partition function for the boundary condition b.c.2 can be expressed as follows: 



1(3 r°={Tf .«..Tf}CA| rf={ 7 f,...,7|}CA| 

fle-MMW x f[ e -^ u U\ } ( 3 .2i; 

j=l k=l 



where 



E 2 (IUT) := J2 ®b(IUT)+ £ $ B (/ur)- £ $ B (/U0) 

B:|B|>3 B:|B|>3 B:|B|>3 

snr^0 Bn/Vd 

= £ $ B (/ur)+ X! (*B(/ur)-$ fl (Ju0)) 

B:|B|>3 B:|B|>3 

Bnr^0 snr^0 
snf=0 

= £ $ B (iur)+ J] %( J "ur), (3.22) 

B:|B|>3 B:|B|>3 

Bnr^0 Bnr^0 
Bn/=0 

where we have defined the quantity 

$' B (J U T) := $ B (J U T) - $ B (J U 0). (3.23) 



It satisfies the bound 



$s(/UT)| < 2c 2 (^r^\ (3.24) 



m(B) 



where m(B) is defined through ( |2.37f ). 



We prove in Sections § and |5] that the 100 and 111 interfaces are rigid at low temperatures. 
The term dependent on the area of the interface, in the expression ( |3.5| ) for the "bare" energy 
£?i are (J) of the interface, is responsible for the rigidity of the 100 interface. However, for 
the 111 interface, the corresponding area-dependent term in the expression for the "bare" 
energy i?2 are (-0 of the interface is not sufficient for stabilizing it against thermal fluctuations. 
The rigidity at low temperatures results instead from the geometry-dependent contribution 
of the plaquette potential h p (I U 0) to the energy of the interface (i.e., the third term on the 
RHS of (CT)). 

The proofs of these results involve an analysis of the convergence properties of the parti- 
tion functions S^ om , S^ 00 and S^ 11 in the limit A /* Z 3 . A direct application of the method 
of cluster expansion requires the contours to be non-interacting. This means that the energy 
of a configuration is given by a sum of terms, each depending on only one contour in the 
corresponding compatible family. This is not true for the FK model, because the long range 
interactions in its effective Hamiltonian induce interactions among the contours. We over- 
come this technical difficulty by rewriting the partition functions in terms of configurations 
of non-interacting but more complicated contours called decorated contours. In the next 
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section we define the decorated contours and derive expressions for the partition functions 
in terms of them. 

3.1 Decorated contours 

Let us first explain how we express the partition function of a system under homoge- 
neous boundary conditions fl2.19|) in terms of decorated contours. The partition function for 
homogeneous b.c. is given by (|3.2j) , which we repeat here for convenience. 



^hom 
"A 



e 

r={7i,...,7n}CA i=l 



-PJi{U)\li\ 



n 

B|>2 

flnr#0 



-j3*s(T) 



(3.25) 



By convention we treat empty products as unity. The idea is to analyze the effect of the 
dominant part Hq 2 ^(U) of the relative Hamiltonian ( |2.23p by a low temperature expansion 
in terms of its contours 7, while treating the contribution of the long range tail, Rx , by a 
high temperature expansion [ |3"T| , |T2"|] . Hence, we write 



^hom 
"A 



rcA 7 er 

Elle 

rcA 7 er 



-0MU)h 



il n [( 



D + i 



snr^C 

|B|>3 



i+E e n( 



(3.26) 



n>l 



s i nr 7 ^0 

|>2 ; (i=l...n) 



1=1 



To each term on the RHS of ( |3.26j ) we can associate a finite set of compatible contours and 
a set of bonds, each bond intersecting the support of at least one contour in the set. More 
precisely, a decorated contour is defined by the pair 



D = (T D ,B 



D 



(3.27) 



where C T is a finite set of compatible contours and Ed := {Bi, . . . B n } is a finite set of 
bonds such that for each Bi, 1 < i < n, there is a 7 G with Bi H 7 7^ and 

(3.28) 



suppD := (U 7eri3 7) U (UBeBr,- 8 



is a connected set. In ( 3.28 ) suppD denotes the support of a decorated contour D. In the 
sequel we shall use the symbol D for both a decorated contour and its support. Moreover, 



\d\ ■= E N + E 9(B). 

76D BeD 



(3.29) 



Decorated contours have the following properties: 



The interiors of any two distinct closed Ising contours, which belong to a decorated 
contour D, do not intersect. 
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• Each bond B in a decorated contour intersects at least one contour in D. 

• Any two contours in a decorated contour D are connected through bonds and other 
contours in D, i.e., for each pair of Ising contours 7, 7' in D, there are contours 

71, 7fc in D such that 7 ~ 71, 71 ~ 72, • • • , 7fc ~ 7', with the understanding that 

two contours 71 and 72 are connected, denoted by 71 ~ 72, if one of the following holds: 

— there is a bond B in D such that B intersects both 71 and 72 

— there are two bonds B\ and B 2 in D such that Bid B 2 ^ 0, Bi intersects 71 and 
B 2 intersects 72. 

To each decorated contour D = {Y Di Bd) we can associate a weight W(D) as follows: 

W(D) := [] e-^WM ]J ( e ~^ B " !)• ( 3 - 30 ) 

7G-D ssd 
|B|>3 

Let T> denote a finite family of compatible decorated contours, i.e., a finite set of mutually 
non-intersecting contours. Then the partition function can be expressed as follows. 

sjr = e n w(d). (3.31) 

where \W(D)\ < W {D) for all D G £>, with 

WoP):= II e-^^-^l J] [exp(/?c 2 (^)^)-l]. (3.32) 



7Gr D Bee D 



The above bound follows from the Peierls bound ( |2.30|) and the exponential decay (|2.32|) . In 



( P-31| ) the partition function for the lattice model under homogeneous boundary conditions 
has been expressed in terms of a gas of non-interacting, pairwise disjoint, decorated contours. 
The methods of cluster expansion can now be applied to analyze its convergence properties. 
The partition functions corresponding to mixed boundary conditions fl2.20|) - ( |2.21| ) can 



be expressed in terms of decorated contours in a similar manner. However, under these 
boundary conditions there is a contour - the interface - which is pinned to the boundary 
of the volume A. In our definition of decorated contours the interface is treated differently 
from the remaining Ising contours in the volume A. It is not considered to be a part of a 
decorated contour. Consequently, in the expressions for the partition functions under mixed 
boundary conditions, there is an additional sum over all possible interfaces. As described 
in Section ^ an interface I divides the volume A into two subvolumes A" and A*. For each 
interface, the set of non-interacting decorated contours, corresponding to a mixed boundary 
condition, can be decomposed into three subfamilies which consist, respectively, of contours 
which intersect /, which are above I and which are below /. Hence, for a given interface 
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/, under the boundary condition b.c.l ( |2.20| ), we define the following compatible families of 
decorated contours. 





■= {{D h . 


Vj 


■= {{D mi 


V] 


■■= {{D ni 



D Tlr \D rH n TV 0, D rH n / = 0, 1 < i < r, rn E Z}}, (3.33) 

where T a and T b are the subfamilies of compatible (Ising) contours which lie entirely in the 
subvolumes and Aj respectively. For a given interface /, under the boundary condition 
b.c.2 ( |2.21| ), the corresponding families of compatible decorated contours are denoted by T>j, 
Vf and V\. 

From ( p7T3|) , Q) and ( gig ) it follows that 



"A 



lei 



e n w(d))( e n 



(3.34) 



where is defined through (|3~30|) and satisfies the bound \W(D)\ < W (D), with W (D) 

being defined through ( |3.32| ). 

For DgP/ 

:= [] e~^ u)M J] ( e "^ - !)• ( 3 - 35 ) 

7G-D 



BSD 
\B\>3 



From the Peierls bound ( |2.30| ) and the estimate ( |3.11| ) it follows that | W/(Z?)| < Wf where 



Wj(D) := J] e-to^W ^ [ exp (2/3c 2 (|)^) - 1 



BSD 
|B|>3 



(3.36) 



Similarly the partition function S^ 11 , defined through (|3.21 ) - (|3.22|) , can be written as 

n w i( d )] 



^111 

"A 



E' 

/ex 



X 



e n w (d)){ e n «m 



D|nA^0 Devi 



D°nA^0 Dex>° 



B-l). 



where for D e T?j 

W f (D) := J] e -^WlT*l J] ( e 

7GD sen 

|B|>3 

The bounds (|2T30|) and ( ggp imply that |W>(D)| < Wj? where 

:= II e^ CoC/ " 1171 II [^MWc2Cjf) m{B) ) 

yeD B&D U 

|B|>3 



(3.37) 



(3.38) 



(3.39) 
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where m(B) is defined through (|2.37| ). 

In order to proceed we need to analyze the convergence properties of series of the form 

Sa:= E II W(D) (3.40) 

DnA^0 Dev 

where T> is a finite set of compatible decorated contours with weights given by 

W(D) := [] e~ m II ( e_/3GB - !)» ( 3 - 41 ) 

7G-D Be-D 

where the function Gb is given by ( p. 24! ) for the homogeneous boundary conditions (|2.19| ) 
and by the functions $' B ( |3.8|) and $' B (|3.23|) for the mixed boundary conditions ( |2.20| ) and 
( |2.21| ) respectively. The bounds satisfied by these functions are respectively given by ( |2.32| ), 
( |3.11| ) and (|3.24|) . Since the series in (|3.40|) is expressed as a sum over compatible families of 
non-interacting decorated contours, its convergence properties can be studied by the method 
of cluster expansions. The convergence of the above series (for Gb = $b and Gb = &b) 
follows from Lemma [| given below. 

Lemma 5 Consider the series 

Sa := E II W(D), (3.42) 
x>nA^0 Dev 

where 

W(D) := [] z~ mi) II ( e_/3GB - !)> ( 3 - 43 ) 

7GD BeD 

and assume that there exists positive constants C\ and C2 such that 

E(i) > CiA| 7 | (3.44) 

and 

\G B \<C 2 \^ B \ (3.45) 



where g(B) is defined through ( 2.1C ) and A < 1. Then there exists constants bo, Xq > such 
that for f3\ > bo and A < Ao, the series has a convergent cluster expansion, i.e., 

logS A =E^ E ■•• E * T ({ Dx, - - - ,D N }) (3.46) 

N>1 JV - DinA^0 D N nA^9 

(the cluster expansion,) ; where \1/ T is a function on families of decorated contours with the 
property that 

y T ({D u - ■ ■ ,D N }) = if {D ir --,D N } is not a cluster, (3.47) 
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i.e., if D\ U • • • U -Djv is n °t a connected set. It satisfies the bound 

\^(D U ...,D N )\ 

^ UTTi — mTT - SN (3 - 48) 
{ Di ,...,d n }30 l^i u • • • u D m 

where Sn is a constant of the order of 

N 

supniwpoi, 

i=l 

i/ie supremum being taken over all sets of N decorated contours (not necessarily pairwise 
disjoint). 

The above lemma is proved in Appendix B. The proof of the convergence of the series S\ 
( P-40| ) for Gb = &b is analogous to the proof of the above lemma and is hence not included. 
The only difference is the replacement of the bound ( 3.45Q by the bound 

\G B \ < C 2 \ m[B) . (3.49) 

where m(B) := max(5, g(B)). 

Let P denote a cluster of decorated contours. It is a connected set of intersecting deco- 
rated contours. A single decorated contour can occur several times in a cluster. Further we 
define 

1^1 ==£ |£>l =£ (EM +£</(£)) (3-50) 

DeP DeP 7G-D BeD 

Let \I/ T (P), tyJ(P) and tyJ(P) denote truncated functions, defined on the cluster P, which 
satisfy the following bounds. 

\^ T (P)\ „ 
E^rn ^ sup]l\W(D)\, (3.51) 

< suplli^p)!, (3.52) 

and 

^ su P niW>(I>)|. (3.53) 

If the cluster P consists of N decorated contours then the product is over a set of N decorated 
contours and the supremum in the above estimates is taken over all such of N decorated 
contours. 

To determine whether the 100 and 111 interfaces are rigid, we need to analyze the fol- 
lowing quotients. 

^100 ^111 

Z 1 J°° = ^ f - and Zi n = =^. (3.54) 

"A "A 

Using the results of cluster expansions (see, e.g., |3H|| ) we can express these quantities in 
terms of the truncated functions defined above. 
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Proposition 6 There exist three constants Uq, B\, and B 2 , independent of the volume A, 
such that for all U > Uq and and (3/U > B\, the quotient Z\ 00 can be written as follows: 



rvlOO 



E 

/ex 



exp £ *f(P)- £ * T (P) 



P: 
Pn/^0 

PnA^0 



P: 

Pn/^0 
PnA^0 



(3.55) 



For a// U > Uq and (3/U > B 2 , the quotient can be written as follows 



rylll 



/ex 



-/3£j? are (7") 



x exp 



E *?(p)- E 



P: 

Pn/^9 
pnA^0 



P: 

PnA^0 



(3.56) 



Proof: The proof is standard. We expand the partition functions appearing in the numer- 
ator and in the denominator of Z\ 00 and of Z A n by the cluster expansion performed in the 
the Appendix B. The truncated functions which are left are precisely those defined on the 
clusters which intersect the interface. H 
Alternatively Z A n can be written as 



rill 



/ex 



-PE% cc (i) 



(3.57) 



where E^^I) denotes the energy of the interface in the presence of the decorated contours 
and is defined as follows. 



Pn/#0 
PnA^B 



P: 

Pni^B 
PnA^B 



(3.58) 



We shall refer to -R? 60 ^) as the energy of the decorated interface. 

4 Rigidity of the 100 interface 

The study of the 100 interface requires the second order decomposition of the effective 
Hamiltonian of the FK model, in which the second order truncated Hamiltonian H.q (U) is 
the Ising Hamiltonian, with coupling constant J\(U). We point out that the Hamiltonian 
H.q (U) is generated by the second order quantum fluctuations, this means that the rigidity 
of the 100 interface of the FK model is of quantum nature. The proof of the rigidity of 
the 100 interface of the three-dimensional FK model is a generalization of Dobrushin's 
proof of the rigidity of this interface in the three-dimensional Ising model at sufficiently low 
temperatures. 1151 . 
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Proof of Theorem [l|: 

Let ProbA(-0 denote the probability of occurrence of an interface / in the FK model, 
defined on a finite cubic lattice A, under the boundary condition b.c.l ( |2.20| ). It is defined 
as follows: 



P ™bA(/) = exp {-/?£$•»(!)} x exp £ 9j(P) - £ * T (P) (4.1 



Za 



pn/^0 Pni^0 
PnA^0 PnA^B 



To determine whether the 100 interface is rigid, we need to analyze the properties of this 
quantity in the thermodynamic limit. 



We notice that the above expression (|4. 1| ) is similar to the corresponding probability of 



a 100 interface in the Ising model. This is because the leading term in the "bare" energy, 
£?i are (J) (J3.5D , of an interface / is exactly equal to the energy of an interface I for an Ising 
model with coupling constant J\(U). We point out that this similarity results from the fact 
that we have expressed the probability ProbA(J) directly in terms of an actual interface /, 
separating the two coexisting phases, instead of expressing it in terms of a decorated interface. 
(The latter is given by a connected set (/, £>/, Tj), where Bj is a finite set of bonds and Tf is 
a finite set of Ising contours; [see |33]] and references therein].) This leads to a considerable 
simplification in the calculations since an actual interface between two coexisting phases, 
under the boundary condition b.c. 1, reduces to a flat interface orthogonal to n = (0,0, 1) 
at zero temperature. This property is, however, not satisfied by a decorated interface. 

Let us describe the minor differences which arise in the description of the 100 interface 
of the FK model with respect to that of the Ising model. 

• The first difference lies in the contributions of the truncated functions. For the FK 
model these functions are defined on sets of decorated contours, whereas for the Ising 
model they are defined on sets of Ising contours. However, there exists a positive 
constant Uq for such that for all U > Uq, the contribution of the truncated functions 
is exponentially small for the FK model, as is the case for the Ising model, difference 
arises from the fact that the "bare" energy, E^ are (I) [( p.5| )1, of the interface consists of 
terms in addition to the leading Ising-like term Ji(U)\I\. These terms are, however, 
small for for U > Uq, where Uq is a positive constant. Hence they are treated in the 
same way as the truncated functions. 

The proof of the rigidity of the 100 interface of the FK model is similar to Dobrushin's 



proof of the rigidity of the 100 interface of the three-dimensional Ising model ||15|| . The 
proof of the rigidity can be converted into the study of a two-dimensional contour model, 
which resembles an Ising model with long range interactions, in which the ground states 
are the projections of the ceilings on the 100 plane and the contours are the corresponding 
projections of the walls. The rigidity of the 100 interface at low temperatures follows from a 
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Peierls argument on the contours (walls) |Tq| . Taking into account the small modifications 
described above, we deduce that there exists positive constants Uo and D such that, for all 
U > Uq and (3/U > D , the assertions of Theorem 1 (see Section 0) are true. | 

We would like to remark that Theorem [1], stated for the FK model, is in general valid 
for a wide class of lattice Hamiltonians, namely, those which can be expressed as a sum of 
two terms: a dominant nearest neighbour Ising Hamiltonian, and a remainder consisting of 
long-range many-body interactions satisfying exponential decay (|2.12|) . 



5 Rigidity of the 111 interface 

The case of the 111 interface can be treated by using the ideas of Dobrushin but the 
situation is more involved. The main difference between the 111 interface and the 100 



interface is that the second order decomposition of the effective Hamiltonian (|2.27|) does not 
lead to the existence of a unique ground state interface in the 111 direction. The energy of a 
111 interface w.r.t the leading part, Hq 2 \u) } of the relative Hamiltonian, is proportional to 
the area of the 111 interface. Hence, the ground state interface has minimal area. However, 
there are infinitely many such interfaces in the infinite volume limit (Proposition ^ of Section 
|2.3|) . Hence, the ground state of the 111 interface has an infinite degeneracy. Thus, to prove 
the rigidity of the 111 interface we require a more refined decomposition of the effective 
Hamiltonian, namely, the fourth-order decomposition. In this section we prove that the 
degeneracy of the 111 ground state interface is lifted by the fourth-order truncated effective 
Hamiltonian, Hq (U), which takes into account the fourth-order quantum fluctuations. The 
study of the properties of the 111 interface can be reduced to the analysis of a model defined 
on a two-dimensional triangular lattice which is obtained by projecting the interface onto a 



fixed plane (see Sections ^7|] and |5.1| for details). We refer to this model as the rhombus 
model. The rigidity of the 111 interface at sufficiently low temperatures can be deduced 
from the low temperature behaviour of the rhombus model. 

We first define the main quantity required in our proof of the rigidity of the 111 interface: 
the probability of occurrence of an interface J, in the volume A, under the boundary condition 
b.c.2 ' 



Prob A (J) = _L.exp{- / 9E5 are (J)} x exp <J 



e e * T (^) 



P: 

Pn/#0 
PnA^0 



P: 

Pn//0 
pnA^B 



(5.1) 



5.1 Description of the Rhombus configurations 

As we have seen, the projection of each face of the interface /, onto the plane P yields 
a rhombus of one of the three types, i.e., belonging to TZ , 7?i,or 7Z 2 - An interface / is 
projected onto a covering of P with rhombi, which we refer to as an ^-configuration (see 
Section |2.2| ). In general, this rhombus covering is not a tiling. Moreover, many interfaces 
are represented by the same R configuration. Let C\ denote an R configuration on Pa. 
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perpendicular to them. The projection of such a line segment on the plane P is referred to 
as a A-link. A single A-link can be the projection of several such line segments. Hence, each 
A-link has a multiplicity, m, which counts the number of line segments which project onto 
it. 

5.2 The energy of a rhombus configuration 

To order U~ 3 , the bare energy E^^^I) of an interface J, defined by (|3.20|) , is a function 



of the local geometry of the interface, and hence is a function of the ^-configuration. The 
plaquette potential h p , (given by (|3.17|) ), as well as the next nearest neighbour interaction 



h{ x ,z} with \x — z\ = 2, (given by ( [3.18; )), play crucial roles in determining this energy. Their 



contributions to the energy for different spin configurations are given as follows: 

= -16 (5.6) 
= -12 (5.7) 

= (5.8) 

= (5.9) 

&{*,*}(+ - -) = >*{*,*}(- - +) = "2 (5.10) 

The plaquette and next nearest neighbour configurations which correspond to the lowest 
energy for these potentials are those given in (|5.6|) and (|5.10| ) respectively. Moreover, from 




( |5~2"D it follows that the plaquette configuration with lowest energy corresponds to a good 
pair of rhombi. Hence, a connected set of faces of the interface whose projection on the 
plane Pa consists entirely of good pairs of rhombi defines a local ground state configuration 
of the interface with respect to the truncated Hamiltonian [( |2.34|) 1. Since two rhombi 
that form a good pair are necessarily of the same type, it follows that there are three such 
local ground state configurations corresponding to rhombi in the families TZq, 1Z\ and IZ2 
respectively. 

We conclude from the above that those minimal area interfaces whose projections on 
the plane P are tilings with rhombi belonging to a single family have minimum energy 
with respect to ( 2.34j ) and are hence referred to as ground state interfaces. Since we 



have chosen the standard b.c. for the -R-configuration to be given by a tiling of the plane 
P\Pa with rhombi in the family 1Z$ [see Section |2.2|| , it follows that a ground state interface 
projects onto a tiling of Pa with rhombi in the family 1Z$. We hence denote a ground state 
interface as I-ji . It is easy to see that I-ji has a perfect staircase structure (see Figure |5J). 
Thus, under the boundary condition b.c. 2 ( |2.21| ), the geometry-dependent contribution of 
the plaquette potential h p leads to the selection of a unique ground state interface (up to 
translations) from the infinitely many minimal area interfaces. 
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Figure 5: The "perfect staircase" structure of the interface I-n§- 



For any interface / the connected components of the set I\I-r. are referred to as pyramids. 
They represent the local distortions of the interface from a perfect staircase structure. 

We would like to point out that the above mentioned phenomenon of ground state se- 
lection is a consequence of the Fermi statistics of the electrons in the model. The anticom- 
mutation relations for the electron creation and annihilation operators play a crucial role 
in determining the exact expressions for the potentials h p [( |3.17|) 1 and h{ XjZ ] [ (|3.18|) 1 which 



are responsible for lifting the infinite degeneracy. If instead of fermions we consider bosons, 
then the corresponding commutation rules for the creation and annihilation operators yield 
the following expression for the plaquette potential: 

hp 1 S x SyS Z&W ~t~ ^\^x^z ~t~ SySyj 2), (5.11) 

where w,x,y and z are four sites forming a plaquette. It is easy to see that 



> M ~ > ( 5 - 12 ) 





and hence, in the bosonic case, the perfect staircase configuration is not favoured by the 
plaquette potential in fourth order of perturbation theory. 



5.3 The energy of a ground state interface 
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The energy of the ground state interface In is given by 

ErVn ) := £ 2 bare (4 ) - \[ E *l ( p )~ E (5-13) 



P: 

PnA^0 PnA^0 



where, 



^ re (4 ) = (J 2 (f/)- i ^)|4 |-(^ + h.o.)iV 1 (4o)+ E ^(4 O U0). (5.14) 



{B:\B\>3;BnI no ^} 



Here |/ft | is the total number of faces in the interface In , which is equal to the number of 
rhombi in the corresponding tiling of Pa, and Ni(In ) is the total number of shared edges 
of rhombi in the tiling. The above expression ( |5.14j ) follows from ( |3.20|) . It is obvious that 



£Tb are (/-R.J, j s the same for i — 0, 1 and 2. Hence, we can define the relative energy of an 
interface / as follows. 

ef c {I) := Et(I) ~ Et(ino). (5.15) 

5.4 The components of an i?— configuration 

An Reconfiguration can be decomposed into bases - which are local ground state config- 
urations of [( [2.34| )1, and R-contours - which represent the excitations. 



Bases: A base of a given R configuration is a maximally connected set of good pairs of 
rhombi. As each base is the union of open sets, it is an open set, and by definition distinct 
bases do not intersect. Denote by {Ci, .., C p } the family of bases of a given R configuration. 
The type of a base Q is defined as follows: r(Ci) = j, if the rhombi it consists of belong to 
IZj. Due to the mixed boundary conditions one of the bases is connected to the boundary 
and hence becomes infinite in the thermodynamic limit. Let C\ denote this base. Since the 
standard b.c. is chosen to be a tiling of P \ Pa with rhombi in 7?-o, we have r{C\) = 0. 

R-contours: The maximally connected components of the complements of the bases, i.e., of 
PAU^Ci, are called the -R-contours. Isolated vertices are not considered to be i?-contours. 
Each i?-contour is a closed complex and is denoted by the symbol T. 
An i?-contour T is defined by a pair 

T = (suppT,0 T ) (5.16) 

where 

• suppT denotes the geometric support of the -R-contour and is a connected subset of 
the plane P, 

• ©t denotes the configuration on this support. It is defined in terms of overlapping 
rhombi, 5- and u- lines and A-links that span suppT. 
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For notational simplicity we shall often use the symbol T to denote both the i?-contour and 
its geometric support suppT. 

An /^-configuration, under standard boundary conditions, is given by a compatible family 
of i?-contours, Ca := {Y 1 ,...,Y n }, i.e., a finite set of non-intersecting, closed i?-contours, 
along with a specification of the types of the bases separating the i?-contours. To avoid 
complicated notations, the types of the bases will be specified only when required. Those 
i?-contours whose supports are not in the interiors of any other i?-contours in C\ are referred 
to as the external contours of C\. 

5.5 The structure of the R— contours 

Each i?-contour T has a detailed structure. It can be decomposed into two families of 
subcontours, one of which can be empty. 

• The overlapping R-subcontours {T™, Y^} are the maximally connected sets of over- 
lapping rhombi contained in T. Each overlapping i?-subcontour is a closed complex. 

• The complement of the overlapping i?-sub contours in T, i.e., T \ U^lfT^ is consid- 
ered as an open complex, whose elements {Tf* , T**} are maximally connected open 
complexes. They are called the standard R-subcontours. We refer to an i?-contour as 
a standard i?-contour if it has no overlapping components. 

Hence, an R contour T has the following decomposition: 

T := T{p,q) = {UgT?} (J {uftTjf } (5.17) 

Let us describe the detailed structure of the different i?-subcontours. A standard R- 
subcontour T|* is uniquely specified by a configuration A|* of 5-lines. We refer to such 5- 
lines as standard 5-lines. An overlapping i?-subcontour T° v is characterized by four families 
of configurations. 

• A configuration of overlapping triangles T° v . It is easy to see that each overlapping 
triangle in the projection of an interface has an even overlap number and each rhombus 
in T° v contains at least one overlapping triangle. Let \% ov \ denote the sum of the 
overlap numbers of all the overlapping triangles contained in T° v , i.e., 

\m ■= E o(t). (5.18) 

teT° v 

Let a° v denote the number of extra faces of an interface (in comparison with the number 
of faces of a minimal area interface) which project on to the support suppT^. It is 
given by 

\'Tov I 

Since \T° V \ is even, a° v is an integer. 
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• A configuration A°" of 5 lines. The total length of the 5 lines in T° v is denoted by 
\A?\. 

• A configuration tt° v of u lines. The total length of the uo lines in T° v is denoted by 

• A configuration A° v of X links. The total length of the A links of T° v is denoted by 
|AT|. 

The number of rhombi spanning supp T° v is bounded above by af , since each rhombus in 
supp T° v is an overlapping one. Each of these rhombi contribute at most three distinct sites 
to supp T° v since the latter is a connected set. The number of sites in the R contour T(p, q) 
(defined through ( |5.17| )) hence satisfies the bound 

|T| < E{3|af| + (|Af| + 1) + (|AT| + 1) + (|nr I + 1)} + E(|Af | + 1). (5.20) 

i=l k=l 

The relative energy of the interface I can be expressed in terms of the energy of the corre- 
sponding R configuration Ca = {T 1; T n }. It is given by 

i=n 

ef c {I) = ef°{{T 1} TJ) := £ F(T<) + W{{T U TJ), (5.21) 

i=i 

where F(Ti) is the energy of the i?-contour Tj computed with the fourth order truncated 
relative Hamiltonian defined in fl2.34|) , and W({Ti, T n }) is the contribution to the 
energy of the R configuration arising from the long-range tail Rx 5 defined in The 
latter consists of higher order corrections to the energy of each individual contour, as well as 
interaction energies between the different R contours in the compatible family. From (|5.15|) , 
O and Q3758[) it follows that 

n -I 

E*TO = J 2 (C/)(|/|-|4 o |) + ( +h.o.)[ E h {x , z} (I)- £ h M {Ino) 

i=l {i,z}eB? {x,z}e-B? 

1 r K 
+ (T7^ + h -°-) f E MO" E (5-22) 



16t/ 3 

and 



p={i,3/,z,t}ePf p={a;,3/,z,t}ePf 



W({T 1 ,...,T n })= = { £ $b(/U0)- E $b(/*oU0)} 

S:|S|>3 S:|B|>3; 



P: P: 

Pnl#0 Pn/ K 
PnA ^" PnA#0 
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+ \{ E * T (^)- E 



Pn/#0 
PnA/0 



pn;. 



^0 



PnA^0 



(5.23) 



Lemma 7 Let T &e an R-contour defined through ( 5.1% ) and let C\ denote an R-configu- 
ration. The energies F(T) and W(C\), defined through ( |5. 22\ ) and ( \5.23j) respectively, are 
given by the following expressions: 



F(T) = J 2 (f/)gar + ^ 2 (f/)(E>fl+El A r 

i=l 



k=q 



k=l 



t=l 



l=p 



i=i 



4[/ 3 



OtI I 



1=1 



where 



j 2 (u) = -L-iL + h,,.. 



2Z7 8£/ 3 



and 



B:|S|>3 
S p nC A #0 

- ^{ E E 



/3 ,. 

p p nC A ^0 



P: 

PpnC A #0 



where 



(5.24) 

(5.25) 
(5.26) 



(5.27) 



$ B (C A ) := $ B (/ U 0) - $ B (/^ U 0) (5.28) 

and i5p tae projection of the bond B on the plane P. The notation B-p R Ca 7^ is used 
to denote the condition that Bp intersects either of the following in the R configuration: an 
overlapping triangle, a 5-line, an uo-line, or a X-link. Similarly, Pp denotes the projection of 
the cluster P of decorated contours, on the plane P. 



Proof: We use ( |5.6| )-( |5TTU[ ) and the definition of the different components of an P-contour 
to obtain the expression for the energy F(T) from ( 5.22|) . It is a sum of the the energies 
of the overlapping- and standard i?-sub contours in T. The energy of an overlapping P 
subcontour is given in terms of the number of overlapping triangles, 5-lines, w-lines, and 
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A-links which constitute it. The energy of a standard i?-subcontour is proportional to the 
number of 5-lines in it. 

The expression ( |5.27| ) follows from the definition ( |5.23| ), since the only terms which survive 
in each of the three paranthesis on the RHS of (|5.23|) are those in which the projections of the 
bonds B or the clusters P on the plane P, intersect at least one .R-contour in the compatible 
family Ca- This concludes the proof. 

In the definition ( |5.28|) we have made use of the fact that while the interface / corresponds 
to the ^-configuration Ca, containing the -R-contour T, the projection of the ground state 
interface I-ji contains no .R-contours. | 

More generally, the energy of an Reconfiguration Ca can be expressed as follows: 



_dcc 



(C A ) = £ e(T | C A ), (5.29) 
Tec A 



where e(T | Ca) denotes the total energy of an i?-contour T belonging to the Reconfiguration 
Ca- It is defined as follows: 

e(T | C A ) := F(T) + W(T | C A ), (5.30) 
where F(T) is given by ( |5.24j) and W(T \ C\) is given by 
W(T\C A ) = ®b{Ck) 

B:\B\>3 



- ^{ E * t i{p)- E (5.3i; 



P: P: 



5.6 The relevant probabilities 

The probability of occurrence of an interface / in the volume A can be identified with 
the probability of the corresponding compatible family of i?-contours Ca = {T 1; Y n } in 
the rhombus model. 

e -/fed ec ({T 1) ...,T n }) 

Prob A (/) = Prob A ({T 1) T n }) = Zq^r^y) ^ 

Further, let ProbAT denote the probability of occurrence of an i?-contour T. It is given by 

Prob A T = ^ e _ £ ae C(CA) , (5.33) 

where the sum in the numerator is over all Reconfigurations which contain the given R- 
contour T, while the denominator has an unrestricted sum over all Reconfigurations. 

To prove the rigidity of the 111 interface, we need to find an upper bound to the proba- 
bility ProbAT defined in (|5.33|) . This is given in the following proposition . 



Rigidity of FK Interfaces 



39 



Proposition 8 There exists positive constants Uq and bo, such that for all U > Uq and 
f3/U > b , the probability, Prob\(T (p,q)), of occurrence of an R-contour T(p,q) (defined 
through ( \5. 1 TQ ) satisfies the following bound: 

% — t) — q 

Prob A (T(p,g)) < J] {e-^ u ^ v } x {e^ K ^ A ^} (5.34) 

i=l k=l 

where 

HU) = J 2 (U)-A 1 U- 5 -^e 

K 3 (U) = K 2 {\J) - A x \J- h --e-^ MV ~ x (5.35) 

P 

and A\ and A 2 are positive constants depending on C\, c 2 and Uq. The constants J 2 (U) and 
K 2 (U) are defined through (\5.23j) and ( \5. 2dj ) respectively. 

The proof of this proposition requires two steps. The first is to obtain a lower bound on 
the total energy, e(T | C\), of an i?-contour T which belongs to an R configuration C\. This 
energy is defined through (|5.30|) . 

The second step is a generalization of the Peierls argument analogous to Dobrushin's 
treatment of the antiferromagnetic Ising model |14| . A unique specification of an R configu- 



ration requires the specification of not only a compatible family of i?-contours, but also the 
type of the bases adjacent to the inner and outer boundaries of each i?-contour. In other 
words contours in the R configuration are not only required to be pairwise disjoint, but 
there is the additional requirement of matching of boundary conditions. One way of analyz- 
ing such contour expansions would be to use the Pirogov Sinai theory extended to interacting 



contours |J2[ [37||. However, instead of doing this we resort to a much simpler method. We 
use a recipe for removing a contour from a compatible family which is a generalization of 
the one introduced by Dobrushin in the study of the antiferromagnetic Ising model. The 
idea is to map a configuration of i?-contours {Y 1 ,...,Y p } to a new one {T2,...,T'} with 
one contour less, where T r and T' r are either the same, or related to each other by a simple 
geometric transformation. The transformation preserves the energy of the i?-contour, at 
least to order U~ 3 . By using this generalization of Dobrushin's construction we avoid using 
the Pirogov-Sinai theory. 

Step 1: A lower bound to the total energy of an R— contour The following lemma is 
necessary to determine a lower bound to the total energy e(T \ Ca), [defined through ( |5.30|) 1. 



Lemma 9 There exist positive constants Uq, bo, such that for all U > Uq and (3/U > bo the 
following bound is satisfied: 

\W{T\C A )\ < Wo(T|C A ), (5.36) 
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where 



W (T | C A ) := |T| 1 ai U~ 5 + ^e-^^ 1 \ , (5.37) 



with c, a\ and ai being positive constants depending on c\, c<i, Uq and bo- 
Proof: From the definition ( |5.31| ) of W(T \ Ca) it follows that 
\W{T\C A )\ < K(Ca)| 



B:\B\>3 



+ \{ E + E 

|^(Ca) 



P: P: 



< E E 



X *£T B:\B\>3 
B3i 



\B\ 



P3i P3i 



where x* is the projection of the site x of the lattice on the plane P. From the definition 
Q5.281 ) of $ B and the bound (|2.36|) it follows that for all U > c d c 1 , 

E{ E ^f^}<ai|T|£/- 5 , (5.39) 

X*eT B:\B\>3 \^\ 
BSx 



where a\ is a positive constant depending on Uo, c i an d °2- Further, using the bounds (|3.51 
and (|3.53|) , and the definitions (|3.32|) and ( [3.39| ), we obtain the following bound: 

+ ,5,0, 

P3x P5x 

where c and a 2 denote positive constants depending on ci, c 2 , C/ and b . The factor of six 
in the exponent arises from the fact that the smallest Ising contour has six faces. | 
The lower bound to e(T | Ca) is given by the following corollary. 

Corollary 10 There exist positive constants Uq and 6 0; such that for all U > Uq and (3/U > 
bo, the total energy e(T | Ca) f (\5. 3Q) J, of an R-contour T (defined through (\5. lTj )) which 
belongs to an R configuration Ca, satisfies the following bound: 

\e(T | C A )| > MU) J2 < + K 3 (U) £ |Af |. (5.41) 

i=l k=l 
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where 

J 3 (U) = J 2 (U) - Q ai U- 5 - ^e-W^u- 1 (5.42) 

K 3 (U) = K 2 (U)-2 ai U- 5 _ ^e- 6/3a2l/ "\ (5.43) 

and c, a\ and a 2 are the positive constants of Lemma |^. The constants J2{U) and K2(U) 
are defined through (\5.23j ) and respectively. 



Proof: From the definition (|5.30|) of e(T | C\) it follows that 

|e(T|C A )|>Fo(T)-Wo(T|C A ), 

where Fq(T) is a lower bound to F(T), and Wo(T|Ca) is an upper bound to W(T\C\). 
Obtaining F (Y) from (|5.24|) , using Lemma || and the bound (|5.20|) yields the bound (|5.41|) . 



The leading contributions to the energy of the overlapping R subcontours Y° v of the R 
contour T stem from the overlapping triangles T° v of T° v . An edge of an overlapping 
triangle can in general coincide with either one of the following - a 5-edge, an cj-edge, or a 
A-link. A uniform lower bound to the energy e(T | C\) is obtained by omitting the positive 
energies of these additional edges. | 

Step 2: A generalized Dobrushin's transformation. 

Consider an ^-configuration C\ defined by the i?-contours Ti,...,T;. Let Ti be the 
contour we want to remove; = P \ Ti is an open set, which has maximally connected 
components {Oo, O r }, where Oo is its exterior and {Oi, O r } are the components of its 
interior. As Ti is a maximally connected component of the complement of the bases, any 
i?-contour in Oi is not connected to Ti. Therefore, we can uniquely define a type r(Oj) 
of each interior, according to the type of base that separates Ti and the contours in the 
interior of Oi. For simplicity assume that t(Oq) = 0. We would like to lift Ti out of the 
^-configuration and fill the gap thus created with bases of type 0. This can be done only if 
r(Oj) = for every i, which is not the case in general. The extension of Dobrushin's trick 
HU is first to apply a translation Si to Oi such that 

r(S J [O l ])=0,^ = l,...,r. (5.44) 

It is easy to see that such translations Si over one lattice spacing exists; e.g., we can use 
vertical translations over one lattice unit in the upward downward direction. Let us denote 
these translations by S and S~ l , respectively. Then the following relation holds for a base 
Co of type 0: 

T (S n [C }) = -nmod3. (5.45) 

Now we will apply S to each Oi of type 1, and S^ 1 to each Oj of type 2. Let = 0, ±1 be 
the exponent such that: 

r(S n *[O t ]) =0, z = l,...,r. (5.46) 
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In general the translations of the different Oi can now overlap, i.e., 

S^p,] nS n *[Oj] ± 0. (5.47) 

As long as all intersections are bases of type 0, there is no problem, and a new configuration 
with the contour T x removed can be defined. This corresponds to the situation in which 
the i?-contours in S ni [Oi] and S nj [Oj] do not intersect each other. We now prove that this 
is indeed the case. Denote the inner i?-contours in 0, by T^i, . . . , Tj j?i . Then we have the 
following property of the interiors. 

Lemma 11 For every R-contour T it j in the interior Oi we have: 

^[T M ]C0"; j = l,..., qi , (5.48) 

where Oi denotes the closure of Oi, and where rij is an integer such that r(S' ni [O i ]) = 0. 

Proof: As Oi is a simply-connected set, and Tjj is closed, it is enough to show that, 
for every edge e C Tjj, we have S ni (e) C 0«. Let P e be the union of all closed triangles 
intersecting e either with an edge or a vertex. Then we have 

P e C 0~, for all e C T^-, (5.49) 

because otherwise T it j would be connected to Of. It is then obvious that 

S n >[e] C P e C 0~. (5.50) 

■ 

Lemma 12 For any pair of R-contours T il j 1 and T i2 j 2 , the following is true: 

S n ^ [T luJ1 ] n [T i2j2 ] = 0, for all i 1 ± t 2 ; 1 < j x < q h ; 1 < j 2 < q l2 . (5.51) 

Proof: As the i?-contours Ti lt j 1 and T i2j2 are closed complexes in two distinct interiors 
and O i2 , we have 

d^.TO^l. (5.52) 

Hence we have 

d(T iujl: T i2 , j2 ) >2. (5.53) 
As S ni * and S ni * are both translations over one lattice unit, it is clear that 

5"MTn,, 1 ]n^[T l2J2 ] = 0, (5.54) 
whenever ^(Yj^, T i2 j 2 ) > 3. Hence the only case that we need to investigate is 

Let us suppose that 

5"MT njl ]n^[T l2j2 ]^0. (5.55) 
Then there should exist a vertex v such that 

veT,, v e S^[T il:jl ], veS^[T i2j2 \. (5.56) 
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v- 



Figure 6: The relative postions of v, v\ and V2 used in the proof of Lemma 12 



• V 2 



• V l 



• V l 



• V l 



(a) (b) (c) 

Figure 7: The three cases used in the proof of Lemma 



If ri,\ = n 



'2 ' 



then we conclude that Tj lj j 1 and T i2 j 2 have to intersect, which is excluded 



by hypothesis. 



Hence, we must have = 1 and n, 2 = —1 (or = —1 and rii 2 = 1). This means 
that there is v G Ti, V\ G and Vi G Tj 2j2 as Figure ^|. 

The points v, v\ and t>2 are related: S(vi) = v = S' 1 ^). As the three vertices belong 
to non-intersecting i?-contours, there are only the three possibilities shown in Figure 
0. Using the fact that the contour boundaries cannot subtend angles of | when the 
i?-contours do not intersect, we complete these diagrams as shown in Figure || In each 
case we can show that this leads to contradicting assignment of types to the rhombi 
containing the vertex v. For example, in the case (a) we are forced to assign the types 
of the rhombi as shown in Figure |9|. This contradicts the fact that 5* _1 [2] = and 
S(l) = 0, and hence is not allowed. We exclude the case (b) in the same way. For 
case (c), consider the dotted hexagon. It is easy to see that all rhombi with diagonals 
that are edges of the same hexagon must be of the same type. This implies that: 
r(Oj 1 ) = r(Oj 2 ) contradicting the condition n. 



-n 



12 ■ 
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Figure 9: The assignment of types to the rhombi of Figure || (a). 
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Now we can complete the Dobrushin argument and the proof of Proposition |8|. We have 
shown that any i?-contour T can be removed, meaning that the following operations were 
performed to obtain a new configuration of non-overlapping i?-contours. 

• Erase T 

• Translate the interior Oi by S ni such that 

T(S n *[O i ])=T(O Q ) = 0, (5.57) 

where Oo is the exterior of T. After the translations have been performed some parts 
of the bases, which are now all of the type 0, will overlap. 

• Fill up the gaps that were left with the base of type 

The essential point of this construction is that the R contours which lie in the interiors 
of Ti have been translated without intersecting each other. The effect of the translation is 
to modify only the energies corresponding to the potentials of range greater than two, which 
are corrections to the leading terms. This concludes the proof of Proposition || 

5.7 The Peierls Condition for the geometric R— contours. 

There may be many i?-contours which have the same support. Different i?-contours 
whose supports coincide differ from each other in the configuration of their constituent 
overlapping i?-subcontours, e.g. in the overlap numbers of the overlapping triangles, the 
multiplicity of the A-links etc. 

It is convenient to group the i?-contours into equivalence classes depending on their 
support. This allows us to obtain bounds on relative probabilities entirely in terms of the 
supports of the .R-contours. We define equivalent contours as follows. 

Two i?-contours 

T := (USTD U (u£?Tf ) and T' := (U*Z?Tf ) U (u£?Tf ) 
are said to be equivalent iff they fulfill the following conditions: 

• The standard subcontour T|* 6 T is identical to the standard subcontour T' fc s * G T' 
for k = 1, . . . , q. Each standard i?-subcontour is given by a unique configuration of 
5-lines. 

• The overlapping i?-subcontour Tf° G T and the overlapping i?-subcontour T'° v G T' 
have the same support, denoted by supp(Tj ), for i — 1, . . . , p. 

A geometric contour T = T(p, q) is an equivalence class of i?-contours, i.e., 

T(p, q) = T := {T a = (USTS) U( U "=W ) I SU PP T » = su PP T r for 1 < i < p} , 

(5.58) 
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where the subscript a is used to label the different .R-contours which have overlapping 
subcontours of identical support. Like an .R-contour, a geometric contour can also be de- 
composed into overlapping and standard subcontours: 



T(p,g) = (uSlO U (utfT 



=: fuf, (5.59) 

where the symbol T i denotes the i-ih overlapping geometric subcontour. It follows from 
( |5.58| ) that all i?-contours constituting a geometric contour have the same support. 



We can associate a unique number r° v to the support supp(Tj ) which is defined as 
follows: 

r° v := minimum number of (distinct) rhombi in which suppT™ can be decomposed. (5.60) 

For simplicity we shall often use the symbol T to denote both the geometric contour and its 
support. 

Each overlapping i?-subcontour is characterized by a configuration of overlapping 
triangles, each of which is labeled by an even overlap number. The i?-sub contours belonging 
to a given geometric subcontour T i differ from each other in the distribution and overlap 
numbers of their overlapping triangles. Since each of these r° v rhombi is an overlapping 
rhombus, the following bound is satisfied: 

C > rf, (5-61) 

where aj™ is defined through Q5.19Q , the subscript a labelling the different i?-sub contours 
belonging to T { . Let ProbAT denote the probability of occurrence of an .R-contour with 
support suppT, with support through (|5.58|) ) in Pa. To compute this probability we need 



to sum over all possible i?-contours which belong to T. 

Proposition 13 There exist positive constants Uq, d and d\ such that, for all U > Uq, 

(3/U > do and f3/U 3 > d\, the probability ProbA(T(p, q)) satisfies the bound 

Prob A (T(p,g)) < JJe-^-^'T x jf e -^" 3 ^ 1^1, (5.62) 

i=l k=l 

where D\ and Di are positive constants depending on Uq, do and d\. 

Proof: The probability Prob A (T(p, g)) can be expressed in terms of the probability 
ProbA(T(p, q)) that an /^-contour T(p,q) occurs in a given Reconfiguration. Using Propo- 
sition |8] we obtain the following bound. 

Prob A (T(p, q)) = Prob A (T(p, g)) 

{TeT} 
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Next we need to estimate the sum over a in the last line of ( |5.63| ). Each overlapping 
subcontour with a given support is obtained by the projection of a set of faces (not necessarily 
connected) of an interface. The number of faces of an interface which projects onto the 
overlapping i?-subcontour Y° v is equal to {o^+rf"). Since all the overlapping i?-sub contours 
in Tj have the same support, it is possible to find an edge, which belongs either to a base or 
to a standard i?-subcontour, such that it intersects all these overlapping i?-subcontours at a 
fixed vertex. (The choice of such an edge is, however, not unique) . From these considerations 
it follows that we can replace the sum over a by a sum over the variables a^, which take 
integer values and correspond to the distinct sets of faces (each set belonging to an interface 
/) which satisfy the following properties: 

(1) The projection of the (disjoint) union of the faces in each such set is connected to a 
fixed end of a fixed edge in the triangular lattice spanning the plane P. 

(2) Each set has + r° v faces with cij > rf. 

Since the projection of each such set of faces is a connected set of rhombi, at least one 
vertex of each rhombus is shared by another rhombus. Hence each rhombus has at most 
three vertices to which the fixed edge can be attached. The number of ways of placing ai + r° v 
rhombi on the suppT^ , which is connected to a fixed end of a fixed edge in the triangular 
lattice spanning P, is bounded by the number of ways of constructing a connected set which 
consists of a>i + rf" rhombi and is connected to this fixed vertex. This latter number, which 
we denote by N ov (ai), is easily seen to satisfy the bound 

N ov {ai) <m a > +r °\ (5.64) 



with m = 36, by the Konigsberg Bridge Lemma [44]. Hence, 



Prob A (T(p,g)) < £ m x e -0(a i+ ... + a p )MU) x JJ e - P K a <p)x\A«\ 

ai>r° v ,...,a p >r° v k=l 
i=l a;>0 k=l 

It is easy to see that there exists positive constants U , d , such that for all U > U and 
(3/U > do, the geometric series in paranthesis (on RHS of (|5.65|) ) converges. Hence, using 



the definitions (|5^2p and (|Q3D of J 3 (U) and K 3 (U) we find that for all U > U , (3/U > d 
and (3/U 3 > d\ (where d\ is a positive constant) the following bound holds. 

Prob A (T(p,g)) < JJe-^-^'-r x jf e -^" s ^ 1^1, (5.66) 

i=l k=l 



where D\ and D 2 are positive constants depending on Uo, do and d±. 
5.8 Proof of the rigidity of the 111 interface 
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The rigidity of the 111 interface can be expressed in terms of the probability, Prob A (s X0 = 
— l|b.c.2), that a lattice site xq = (xi, x 2 , 23), such that x\ + X2 + £3 > 1/2, is occupied by 
a "— " spin under the boundary condition b.c.2 (|2.21| ). To have s xo = —1, the site x must 
be enclosed either by a pyramid of the interface or by at least one Ising contour. As a result 
we have that 

Prob A ( Sa;o = -l|b.c.2) < Prob A T+ E Prob A7 , (5.67) 



T3x* 



732:0 



where x* Q denotes the projection of the site xq on the plane P A . The first term on the 
RHS of ( |5.67| ) arises when the interface has at least one pyramid whose projection on P A 
encloses the point Xq. The second term arises when the "— " spin at x$ is enclosed by one or 
more Ising contours, whose presence does not lead to a distortion of the interface from its 



perfect staircase structure around the site xq. The above bound ( 5.67 ) is also satisfied by 
the probability Prob A (s io = l|b.c.2^ , for x = (x-y, x 2 , x 3 ) with x\ + x 2 + x 3 < — 1/2. 

Hence, to prove the rigidity of the 111 interface we need to estimate the terms on the 
RHS of ( |5.67| ). An estimate of the first term is given by the following proposition. 



Proposition 14 There exist positive constants Uq and b\ and b 2 , such that for all U > Uq 
(3/U > bi and f3/U 3 > hi, the following estimate is true: 



_E Vroh A T<C e 



-c"/3U- 



(5.68) 



where Co and c" are positive constants depending on Uq, b\ and b 2 . 
Proof: From Proposition [13| it follows that 

E Prob A T < E II e- f>u - lDir ? x J[ e^* ^ 



(5.69) 



T 

r 3x * 



where we have used the following notations: T i denotes an overlapping geometric subcon- 
tour (see ( |5.58 )), while |A|*| denotes the total number of 5-lines in the standard geometric 



subcontour T fc = Yf . Further, an empty product is equal to unity. 



RHS of (5.69) 



v n e ~ pu ' lDir i v x n e -pu- 3 D 2 \A% 



T-r°"uT s 



+ 



-f3U- 3 D 2 \A at \ 



T3x* 



TfeT 



e n 

T:=T "uT S ' 
T°V8 



-pu- 



- 1D ^ V [ n {(e- 



-PU- 3 D 2 \A1 



1) + 1} 



T|*eT 
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+ E 



e 



-/3U- 3 D 2 \A st \ 



j2 n e-^ u ~ iD < v {i+j2 e nc 6 "^ -3 ^ 1 ^ 1 - 1 )} 

T3x* T°"6T°" 



n>l Tf t ,...,T« t fc=l 
T°V0 (i=l,..,n) 

(5.70) 

+ £ e-^-^IA^ (5.71) 

where |A S *| denotes the number of 5-lines in the standard i?-contour T st . 

The term on the RHS of (5.70) is similar to ( |3.2(j| ) of Section |3.1| . This allows us to 



prove the bound ( |5.68| ) by using a method analogous to the one used in that section. Each 
geometric contour T appearing in the sum on the RHS of ( |5.70|) consists of a finite set of 



non-intersecting overlapping geometric subcontours and a finite set of standard geometric 
subcontours such that each standard geometric subcontour intersects at least one overlapping 
geometric subcontour in the set. We can alternatively express each such T as a connected 
set of auxiliary polymers, each polymer being a connected set consisting of only a single 
overlapping geometric subcontour and a finite set of standard geometric subcontours which 
intersect it. To each term on the RHS of ( |5.71| ) we can associate a standard geometric 



contour which is given by a standard i?-contour. These considerations allow us to express 
the RHS of ([5.69 ) in terms of the elements of a more general polymer system, the polymers 



being referred to as spider contours (or S contours, for brevity) and defined as follows: 

An S contour, ( is a finite connected set of geometric i?-subcontours, containing at most 
one overlapping geometric i?-subcontour. Hence, in general 

C:=T™|J{lJSrf}, ( 5. 72) 

where T denotes an overlapping geometric subcontour. 

In particular, an S contour can reduce to a single overlapping i?-contour (with no stan- 
dard part) or to one standard i?-contour. We refer to an S contour which has no overlapping 
part as a standard S contour and denote it by the symbol ( st . The other S contours are 
said to be overlapping and are denoted by the symbol ( ov . Every standard i?-subcontour 
which belongs to an overlapping S contour, (defined through ( |5.72| )), necessarily intersects 
the overlapping geometric subcontour T . 

A given geometric i?-contour T, with support suppT, can be built from a finite family 
of intersecting S contours, For an S contour £ given by ( |5.72|) , we define 

|C|:=r- + ^|Af|, (5.73) 

i=l 
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where r ov denotes the minimal number of rhombi needed to cover the support of T° v (see 
( |5.60| )). For a standard S"-contour ( st we define 



IC 



st I 



IA 



st I 



(5.74) 



which is the number of 5-lines in it. 

We define the weight rj(-) of an overlapping S contour ( ov (defined through (|5.72 )) to be 



n - 



-j3U- 3 D 2 \Al t \ 



(5.75) 



In ( |5.75 ) we use the convention that an empty product is unity. Hence the case in which 
the S contour does not contain any standard part, is included in ( |5.75| ). The corresponding 
weight for a standard S contour ( st = T st is defined as 



From (1577DQ , (|5Tl| ), (|5T75| ) and ( |5T6D it follows that 



(5.76) 



T 

T3! 



EProb(T) < £_ E livid) 

n>l ll - Ci,-;<n _ j=l 

ClU...UCn=T 



< E 



n>l 



n 



where 



(5.77) 



(5.78) 



C := E ^(C)e ICI , 

Ox* 

In obtaining the bound ( |5.77| ), we have made use of Lemma 3.5 of (as in ( |B.8D of the 
Appendix B). The proof of the bound ( |5.68 ) reduces to the proof of the following lemma. 



Lemma 15 For each D' > 0, there exist positive constants b\ and b 2 such that for all 
[3/U > bi and /3/U 3 > b 2 the following bound holds: 

]T ??(C)e lcl < D' (5.79) 

Ox* 

Proof: 

LHSof(5.79)= ]T l(De lrl + E viCV^ 1 (5.80) 

( ov Bx* (°Ox* 
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Let us first evaluate the second term on the RHS of ( |5.80 ). To do this we use the fact that 



the smallest standard i?-contour consists of six 5-lines. Hence, 



E viCV^ 1 <EE e~ PU D2k e k . (5.81) 

C st Sx* k>6 ( st 3x* 

| C si | = fc 

Consider each standard S contour, appearing in the above sum, as a connected graph which 
contains a fixed vertex Xq. The maximum coordination number of each vertex in the graph 
is five, because a 5-line which intersects a given 5-line at a fixed end, can emerge in any 
one of five directions (in the triangular lattice spanning P). Then by the Konigsberg Bridge 
Lemma [33 

RHS of (5.81) < 6 E 5 2h e-? u ~ 3D2k e k . (5.82) 

fc>6 

The factor of six arises from the fact that there can be at most six different 5-lines in the 
rhombus model which contains the fixed site Xq. Let b > be a constant, such that for all 
(3/U 3 > b, the geometric series in ( |5.82|) converges. Then for all (3/U 3 > b, 

E V (( st )e^ < e- 6 ^" 3 (5.83) 

where D4 is a positive constant depending on b. Next we need to evaluate the first term 
on the RHS of ( |5.80|) . This term can be further decomposed into two sums, depending on 
whether the fixed site Xq belongs to the overlapping subcontour of an 5* contour or not. 
Using the definition ( |5.72| ) of an overlapping S contour, we can write 

E v{C)e^= E v(C)e lC ° Vl + E V(C>^. (5.84) 



3 U — OV 



'a 



The second term on the RHS of (|5.84j) corresponds to the situation in which the fixed vertex 



is necessarily contained in a /^-standard subcontour of the overlapping S contour. 
Evaluation of the first term on the RHS of (|5.84| ): 



Let us construct T° v starting from Xq. Once a rhombus which has Xq as one of its vertices 
is chosen, the next rhombus which intersects it can be placed in twenty four different ways; it 
can intersect the first rhombus either at any one of its four vertices or along any one of its four 
edges. Further, for intersection either along an edge, or at a vertex, there are three possible 
orientations of the pair of rhombi. This allows us to consider each rhombus as a vertex of a 
connected graph which contains a fixed vertex. The fixed vertex of the graph corresponds to 
the rhombus which contains the site Xq. There are twelve different rhombi in the triangular 
lattice spanning P which contains a given site. Then, by the Konigsberg Bridge Lemma 



contour Y'™, such that it contains a fixed site Xq, and has r ov = n, satisfies the bound 



4J| the number, iV(Y \ti\Xq), of different ways of constructing an overlapping geometric 

ns a fixed site Xq, 

N(T v \n; x* ) < Yld\ n with d x = 24. (5.85) 
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The overlapping geometric contour T° v has at most (3n + 1) vertices at which a 5-line can 
intersect it. Moreover, from each vertex of T there can emerge at most four standard 
5-lines. Further, each such standard 5-line can correspond to either one of two different 
pairs of rhombi. These considerations yield the bound 

C ov n>\ k>l 

Let b[ and b' be positive constants such that for all /3/U > b[ and (3/U 3 > b' the series in n 
and k converges. Then for such values of (3 and U the RHS of (|5.86|) satisfies the bound 

RHS of (5.86) < e^"^ 6 e"^"^ 5 , (5.87) 

where D 5 and D e are positive constants depending on b' x and b' respectively. 
Evaluation of the second term on the RHS of ( |5.84| ): 



To evaluate this sum we make use of the fact that each standard i?-subcontour in ( ov 
intersects T° v . We construct ( ov starting from the site Xq which now belongs to a standard 
i?-subcontour. If the standard i?-subcontour which contains Xq has m S- lines, then there 
are at most m + 1 vertices at which can intersect it. Moreover, there can be at most 
four standard 5-lines emerging from each of the (Sr™ + 1) vertices of T° v and there are two 
possible orientations of each such pair of S lines. From these considerations we obtain 



E vie 



< ^^v^^m+lj^rflV^^fSx (3n + 1) E 5 2fc eV^ 3D2fc }" . 

m>l n>l k>l 

(5.88) 

Let b x and b" be positive constants such that for all (3/U > b' x and /3/U 3 > b" the series 
on the RHS of ( |5.8S| ) converges. Then, for such values of (3 and U the following bound is 
satisfied: 

RHS of (5.88) < e"^" 3 ^ e"^"^ 7 , (5.89) 
where D 7 and D 8 are positive constants depending on b'[ and b" respectively. Let 

b x := max(6i,6'i), (5.90) 

and 

b 2 := max(6,6',6") (5.91) 
Then from ( ggg) , flBTBTD and (|5"l9|) it follows that for all (3/U > b x and (3/U 3 > b 2 



E r /(C)e ICI < C e"^'^ 3 =: D', (5.92) 
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where Co and c" are positive constants depending on b\ and 62 respectively. Hence D' is a 
positive constant which can be made arbitrarily small by making b\ and bi large enough. 
This concludes the proof. | 

To estimate the second term on the RHS of ( |5.67|) we make use of the result of the cluster 
expansion given in Lemma [f| This yields the following proposition. 

Proposition 16 There exist constants C/1,63 > such that for all U > U\ and (3/U > 63 
the following estimate is true: 

Eftob A ( 7 )<Cie- c > rl , (5.93) 

79:ro 

where C\ } c[ are positive constants depending on U\ and 63. 
Proof: We have that 

E Prob A 7 < E Prob A £> (5.94) 

73X0 D: 

where ProbA-D denotes the probability of occurrence of a decorated contour D in A. 

E v% W(D){E' v IlDievW(D')]} 
RHS of (5.94) < 13X0 = . (5.95) 



J2v YId'&v W{D' 

The symbol Ysd is used to denote a sum over all finite compatible sets V = {D'} of decorated 
contours which are compatible with D, the latter being a decorated contour which encloses 
the site xq. In the denominator we have an unrestricted sum. The result of the cluster 
expansion [Lemma || can be applied to both these sums to yield 

RHS of (5.95) < E W o( D ) e*p{- E4 E * T ({A, • • • , Av})}, (5.96) 

D: N>1 iV • D 1 ,...,D N 

= 97 (<)&(«) 

where (i) & (ii) refer to the following conditions: 

(i) D 1 U . . . U D N = P is a connected set (a cluster), 

(ii) PnD^ 0. 

We have made use of the fact that \W(D)\ < W Q (D), where W (D) is given by (jQ2p . 

RHS of (5.96) < E W (D)e W {-\D\ E ^E^rP} (5-97) 

o: V>l iv - P r I 

Ml — 

where 

1^1 = E(l7l+ E 

7SD BO7^0 
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and the symbol X)p denotes a sum over all clusters P which contain a fixed site. The bound 
( |5.93| ) is then obtained by making use of ( |3.32| ) and the bound (|3.51|) on the RHS of (|5.97|) . 



From the Propositions |14| and [L6| it follows that for xq = (xx,X2, X3) such that X\ + X2 + £3 > 
1/2, we have the following upper bound on the probability ProbA(s X0 = —1): 

Lemma 17 There exist positive constants Uq,Do and D' such that for all U > Uq, (3/U > 
Dq and (3/U 3 > D' the following bound is satisfied: 

Prob A (s, = -1) < {C e- c " p/u3 + C ie - c "P /u } forx = (x h x 2 , x 3 ) andxi + x 2 + x 3 > 1/2, 

(5.98) 



5.9 Proof of Theorem g 

We now have all the estimates necessary to prove our main result on the rigidity of the 
111 interface, i.e., Theorem ^. 

Proof: For x = (xi,X2, X3) with Xi + X2 + X3 > 1/2, we have that 



< s x >[b. c .2]= 1 - 2 x lim Prob A (s :c = -l|b.c.2) 



(5.99) 



where x* denotes the projection of the site x on the plane P. 

Lemma [17| provides an upper bound to the probability ProbA(s x = — l|b.c.2) which is 
uniform in the volume A. Introducing this bound, ( |5.98 ), in ( |5.99| ) yields (|2.40|) . 

Similarly for x = (xi,X2, £3) with X\ + x 2 + £3 < —1/2, we have that 



< s x >[b.c.2]= — 1 + 2 x lim ProbA^ = l|b.c.2 



(5.100) 



which reduces to ( 2.41| ) by the analogue of Lemma 17 for ProbA^ = l|b.c.2J with x\ + 
x 2 + x 3 < -1/2. I 



Appendix A: Bound on the remainder term of the ef- 
fective Hamiltonian 

In this appendix we present a complete proof of a bound on the remainder term for the 
effective Hamiltonians derived from the circuit representation of |33|, which is missing in 



this reference. A proof of a similar bound for a more general class of Hamiltonians will be 
given in [[32" |. The bound is essential to control the temperature dependence of the effective 
Hamiltonians. Therefore, we present the proof here in considerable detail, although we 
certainly do not claim that the proof is new: the proof of Lemma ^ (given below) closely 
follows ideas of ||, and, in general, our discussion follows the lines of |33j and fill . 
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A.l Definitions 

The Hamiltonian defined on a finite lattice A C Z d , d > 2 is given by 

H A = H + tV - ri e N e - mNi, (A.l) 

where 

H = 2Uj2W(x)clc x (A.2) 
x-eA 

and 

V = ~J2 ctc y + h.c, (A.3) 

<xy> 

where <xy> denotes a pair of nearest neighbour sites on the lattice. Let us restrict our 
attention to the neutral case at half- filling, i.e., \i e — — U. In order to make the origin of 
various terms in the series expansions more transparent, we do not set t — 1 in this appendix, 
unlike in the main text. 

The effective hamiltonian (for a given configuration S = {s x } of ions) is defined through 
the relation 

exp[-/3# eff (/3, S)} = Tr^ e exp[-(3H A ] (A.4) 

where the trace is over the electronic Fock space (denoted by jF e ). 
Iterating Duhamel's formula 

e -ft Ho +tV) = e -/3H + ft dre - {0 - T )H o (_ tV j e -r(H +tV) ? ^ 

Jo 

we obtain the Dyson series 

n JO JO JO 



e 

n=0 ' 

e -(T n -Tn-l)H (_ t yj . .. e -fr-*i)H (_ tV ^ e -nH (A.6) 
Hence 

e -PH A . = e / 3C /(7V e+ 7v I )g t n £ f dr n T dr n ^-- - T dr x 



n=0 <Xiyi>,...,<XnHn> 

e (^r n )H j A . . . e-^.^^. (A.7) 



As in [|33| we introduce the set C(A) of (classical) configurations associated to the electron 
subsystem. An element X G C(A) is a finite sequence X = (xi, . . . ,x m ) of distinct sites in 
A. The state \X) G !F e , associated to X G C(A) is defined as follows: 

|X>:=[ C t(x 1 )...ct(x 1 )] d |0>, (A.8) 

where |0 >G C(A) denotes the vacuum. 
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The symboH denotes a total ordering of the sites in A, chosen to avoid ambiguities in 
the definition of the phase in (|A.8 ). For convenience we choose the spiral order |J for d = 2 
and an analogous ordering for d > 3. This ordering is chosen to have the property that, for 
any finite set X C A, the set {X'} = {iG x ^ X} of lattice sites which are smaller than 
X, or belong to X, is finite. 

For a given sequence of pairs of nearest neighbour sites, (<Xi?/i >,..., <x n y n >), let 

€j\Xj >= 4 c v , > for 1 < j < n, (A. 9) 

where each Sj = ±1. [Note that \Xj\ = \Xj~x\ for j = 1, . . . , n.} 
Setting r = 0, r n+ i = /3 and X = X n = X we obtain 

oo n 

Tr^e"^ := e pu ^ +N ^ Y.t n {\{{X j \e^ T ^- T ^ Ho \X j )} 

n=0 j'=0 

n 

x (nli^UX) (A.10) 

j=0 

Note that 

H \Xj >= e(X j )|X j >, (A.11) 

where 

e(Xj) := 2C7 x [number of sites in Xj for which W(x) = 1, (i.e., s x = 1)] (A. 12) 

As in |3^| , we introduce the notion of trajectories. Let 7$ be a positive integer variable (with 
r = and r n+ i = /3) which we refer to as the "time". A trajectory ( = C( T i) is a sequence 
x(tq), x(ti), . . . , x(r n+ i) of sites in A such that either 

n ^ r i+1 with x{n) = x(r i+1 ) 

or 

Tj = r «+i with ar(ri) and x(r i+1 ) being nearest neighbour sites on the lattice. 
This last case we describe as a jump. Let 

J(() := {(x(Ti) x(r i+1 )) \n = T i+ i] x(Ti), x(r i+1 ) G C} (A.13) 

denote the set of jumps in the trajectory (. Let T = {(} denote a set of non-intersecting 
trajectories. If g is a function on the trajectories, then we define a "sum" over sets of 
trajectories as follows: 



4 



9(0 ■= 1 + E E E I[(i),(ii)l 



C n>l (B lr ..,R ) (x ,...,x n ) 

C, V_l -n/ Xl £C(A) 

x / dr n ■■• r dT ig (Q . (A. 14) 

JO JO 

where = {(xj denotes a set of nearest neighbour sites on the lattice. By l[E] we mean 
the indicator function of the event E; in particular, l[(i), (ii)] in ( A.14 ) vanishes unless 
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1. the relation (|A.9| ) holds 

2. Xo = X n . 

Hence the RHS of ( |A.1U|) can be expressed as a "sum" over sets of non-intersecting trajecto- 



ries: 

exp[-/3# eff (/3,S)] = e^jfe^e(T) ]J e^ ( n 0. ( A -!5) 

where |T| is the number of trajectories in T, and |£| is defined as follows: 

|C| := (vertical length of C) H ({x : s a . = +1} x [0,/3]), (A.16) 

and 

s(T) := (X|T I] ( II 4(. i+1 l*> = ±1- (A.17) 

CeT (a ! (Ti)a:(T i+1 )>eJ'(C) 

The symbol T denotes that the product is "time-ordered", and e(T) denotes the sign of the 
permutation of the electrons under the action of T. 

A. 2 Circuit Representation 

To a given ion configuration S and a given set of trajectories T we associate a set, Q, of 
oriented circuits as follows: 

Vertical segments of trajectories located on sites x G A with s x = +1, will be considered 
as up-oriented components of circuits; vertical segments of the complement of the set of 
trajectories in A x [O,0\, located on sites with s x = —1, are considered as down-oriented 
components of circuits. On each horizontal bond at which a jump takes place, we draw a 
segment with an arrow in the direction of the jump. The vertical segments together with 
the horizontal jumps form oriented closed circuits. More precisely, an oriented closed circuit, 
u, is a maximally connected component of the oriented segments of the trajectories. Let 
Q = {ujx . . . u n } denote a finite set of such circuits. The space of all circuits compatible with 
an ion configuration S is denoted by W(S) [compatibility means that segments of circuits are 
oriented upwards if they are located on sites occupied by ions and are oriented downwards 
otherwise] . 

Let |a; | be the total length of the vertical segments of a circuit u>, and let p(u) E Z be the 
winding number of to. The latter can take positive or negative integer values, with negative 
values indicating that to winds around the "time"-axis [0,/3] with a downward orientation. 

The following two relations hold: 

EKI = |Z(M+/3pM) (A.18) 
reT z wen 

|T| = |A|-JVi+ (A.19) 
wen 
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As we are concerned with the half-filled case, we always have that X^eoPO^) = 0. 
From ( [AloD , (|A7T8D and (|A~T9|) we have that 



e -i3H eS (p,s) _ e pu\A\ 



: n *), 



(A.20) 



where denotes a "sum" over all sets, Q, of non-intersecting closed circuits (defined 
analogous to QA.14Q ) and e(Q) is the sign of permutation of the electrons under the action 
of the circuits. The following lemma is crucial for the next step. Its proof closely follows M, 



and ideas from 46 



Lemma 18 There exists a function e(u) such that VTL' 



So we can write 



e -l3H cS (/3,S) = e PU\A\ 



jf n *(«o 



with the weight, z(u), of a circuit, to, defined as follows: 

z(uj) =e(uj)e- UH t^\ 



(A.21) 



(A.22) 



where j denotes the number of jumps in u. 

Proof of Lemma [TB] : Let n x (r n ) denote the number of electrons at the site x at the time 
r n . Then the set of sites x G A for which 

W(x)+n x (r n ) ^ 1 (A.23) 

is said to define the defect set, D n (S), for a given ion configuration S, at the time r n . 
The section, T n (cj), of a circuit cj, at a time r n , is defined as follows: 



{x G [u fl (the plane 



r = r, 



.)]}■ 



(A.24) 



In particular, T (u) is referred to as the initial ("time-zero") section of the circuit u, since 
t = 0. It is clear that each x G r n (u;) belongs to the defect set D n (S) (for n > 0). 
Let |r n (u;) > denote the state in the electron Fock space, T e , defined as follows: 



|r»>:=( n 4)x|0>, 



(A.25) 



i6r„( u ) 

s.t. s x = l 



where |0 > denotes the vacuum. More generally, for a set, fl, of circuits, we define the 
section at time r n to be 



r n (Q) := {x G [Q n (the plane r = r n )]} 



(A.26) 
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and the corresponding state to be 

|r B (n)>:=( I] 4)J0>, (A.27) 

s.t. s x —l 

The following relations hold: 

(a) 

\Fi(u) >= $b. |ri_i(w) > for 1 < i < n, (A.28) 

where I?j = {< o^t/i >} denotes a finite set of pairs of nearest neighbour sites, and $b. 
denotes the corresponding set of operators 

:=R%}. (A.29) 

(b) 

r M = r n (co) (A.30) 

The above relation is due to the periodicity in the "time" direction. 
Each circuit u is uniquely determined by the following 

• its initial section T (u) (and hence on the state |r n (u;) >); 

• a sequence of operators $s . , . . . , &s n ', 

• a sequence of "times" r 1; . . . , r n at which these operators act. 
The sign e(u) of a circuit u> is given by 

n 

s{uj) = J] < r <HI $ B i l r i-iH > 
i=i 

= <r H|$B n ---$B 1 |r H > (A.31) 

The second line follows from (a) and (b) above. 

To prove the factorization property it is enough to consider a set, Q, of circuits (compat- 
ible with an ion configuration S) that can be divided into two mutually non-interesecting 
("time" -periodical) subsets, denoted by ujb and ujc- Each subset can be made up of several 
circuits. To each of f2, uo B and uo c is associated a sign. 

The sign of oj b is given by 

e{uo B ) =<T {u B )\&B n ---$B l \To(u)B) > (A.32) 

The sign, e(uc), of the component ujq is defined in an analogous fashion, with the 
sequence of operators • ■ ■ replaced by $c • • • . The set Q is defined by an initial 
section F d(= To(^d)) and a sequence of operators ^d x , ■ ■ ■ , &D n m , which is a uniquely 
determined permutation of the sequence . . . , $s n , • • • , ^c: m - Its sign e(Q) is also 
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defined by ( |A.32| ) with the obvious changes. We need to show that the sign for Q factorizes, 
i.e., 

e(Q) = e(oo B )e(uj c ) . (A.33) 

By iteration of the argument we obtain the factorization into the signs of the individual 
circuits. 

Let us first discuss the computation of the sign of a single component u B . In terms of 
the vacuum |0), the expression of the sign e{u B ) takes a very simple form. To obtain it, we 
observe that 

(r 0B i$| B ---$| i |r B> = (o|c I ^*j B ...$| i cf ofl |o> > (A.34) 

where Cf QB is a product of creation operators that creates the section r os (= T (uj b )), i.e., 



Cro fl :=( II 4).- (A-35) 

*6 r 0B 
s.t. s x — l 

The product 

A B := Cr os $B n ■ ■ ■ * a C * ToB (A.36) 

is a monomial in creation and destruction operators that can be combined into number 
operators, because of the required periodicity, 

A B |0> = e(u B ) |0) , (A.37) 

of u B - The factorization of signs is a consequence of the fact that this combination can be 
made in a well-defined fashion which is not affected by the presence of other (compatible) 
circuits. 

We choose the following procedure which we refer to as circuit collapsing. Consider 
the string of operators appearing in the product A B . For brevity we shall refer to each 
appearance of a destruction or creation operator supported on a site x as an occurrence of 
x. We denote by S(A B ) ("shadow" of A B ) the set of sites x occurring in A B . We start with 
the leftmost operator in the product A B ( |A.36| ). In order to yield a non-zero contribution 
this has to be a destruction operator with support on some site x, i.e., c x . We now move 
this operator through the operators present to its right (i.e., downwards in time), using 
the anticommutation relations, until we encounter the next occurrence of the site x in the 
product ( |A.36| ). This is necessarily a creation operator, cj., since otherwise the successive 



actions of these two operators would yield zero. Hence we obtain a factor 

c x cl = l-n x (A.38) 

times a phase a x lS) (u! B ) (n x being the number operator for electrons at the site x). This 
phase arises due to the anticommutation of the initial c x with intermediate operators. All 
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the operators which appear to the left of this factor are not supported in x and, hence, we 
can move this factor to the leftmost end of the string to obtain 

A B = o£\w B ) (1 - n,) A B , (A.39) 

where A B satisfies 

A fl |0) = a x l \to B ) e{u B )\U) , (A.40) 

and has two fewer occurrences of x. 

Next, we repeat the above procedure for the string of operators defining A B and continue 
pulling out, in the same way, successive phases and factors 1 — n x . Once all occurrences of 
x have been dealt with, we obtain a product of factors (1 — n x ) k = 1 — n x and an overall 
phase a x (u B ), so that 

A B = a x {u B ) (1 - n x ) A B , (A.41) 

where A B satisfies 

A B |0) = a x (LU B )e(uj B )\0) , (A.42) 

and has no occurrence of the site x: S(A B ) — S(A B ) \ {x}. 

We then repeat the whole procedure for A B . At the end, once all the sites x in S(A B ) 
have been exhausted, one obtains that the whole of A B has "collapsed" into factors 1 — n x 
times a numerical factor. A simple replacement in ( [A. 37] ) shows that this factor must equal 
e(uj B ). That is, 

A B = e(u B ) Hi 1 - n *) • ( A - 43 ) 

We are now ready to prove the factorization property that we need: Let Q be a family 
of circuits and let it be decomposed into two ("time-periodical") subfamilies of circuits, u B 
and io c- Then 

e(n) = s(uj b )6(u c ) . (AAA) 
We can write u>c via an operator 

A c := C roc $c m • • • ^ Q oc . (A.45) 

[with T oc = T (ujc)], such that 

A c |0) = e(u c ) |0) . (A.46) 

We shall use the following two consequences of the compatibility (i.e., mutual non- 
intersection) of the two components uj b and u c : 

(CI) The monomials Cp os and Cp oc have disjoint support, so that 

\T 0D ) = 8C* rn C* |0), (A.47) 
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where 5 is a phase factor. Hence, 

e(Q) = (0\A D \0) (A.48) 

with 

A D := C roc C ToB ^n n+m ■ ■ ■ $ Rl C? ofl C* Tqc . (A.49) 

(C2) The occurrence of a creation operator c x in factors of Ag implies that the site x 
becomes, or continues to be, part of the support of ujb at least until there is a further 
occurrence of a destruction operator c x . In particular since ojb and ujc do not intersect, 
we have the following property: 

Between an occurrence ofc x in factors of As and the preced- 
ing (ie. immediately to the right) occurrence of c x in factors ^ 
of Ab, there cannot be an occurrence of x in factors of Ac- 



Moreover, as the whole set Vt is periodic in the "time" -direction, we have that each occurrence 
of c x in Ab must be preceded by an occurrence of c' x . Combining this observation with (C2) 
we obtain the last property needed: 

(C3) Between an occurrence of c x in factors of Ab and the immediately preceding occurrence 
of c x in factors of Ab, there is an even number of occurrences of x in factors of Ac. Of 
course, these occurrences correspond to alternating creations and destructions. The 
same property holds after the last occurrence of c x in factors of A# and before the 
first occurrence of c x in factors of Ag. 

From (CT)-(C3) we conclude that the "collapse" of A B gives exactly the same factor 
as in the absence of the component uq- Indeed, the last occurrence of a site x in factors of 
As is a destruction operator, that we can displace up to the previous occurrence to produce 
a factor 1 — n x . Between these two occurrences there is no occurrence of x in factors of Ac 
because of (C2). Hence c x commutes with all the operators $c encountered during the 
displacement (recall that such operators are monomials of even degree). Thus, the phase 
acquired during this displacement only depends on the operators in A^ and hence it is the 
same phase oJ£\ujb) obtained when collapsing the component lob in the absence of any other 
circuit. Moreover, by (C3) the operators in Ao located to the left of the factor 1 — n x , 
obtained in the above manner, involve an even number of creation and destruction operators 
supported in x. Therefore, we can freely move this factor 1 — n^, all the way to the left to 
obtain 

A D = a£\u B ) (1 - n,) Ad , (A.51) 

where A D has two fewer occurrences of £ in factors of A B but otherwise satisfies (C1)-(C3). 
Iterating this process we collapse A B exactly as done in (|A.39 )-( |A.43 ). We obtain 

A D = e(u B ) II (1-n^Ac (A.52) 
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Combining this expression with ( A.48|) and ([A. 46 ) we get the desired factorization ([A. 44 ). 



A. 3 Cluster expansion 

The logarithm of 



4 n *m 



can be developed in a cluster expansion. Adapting Theorem 3.1 of [36] to our case (in which 
one variable, the "time", is continuous) we have the following result: 



if 4 



\Z\UJ 



|e |w| x[(x, r) G u] < Cj with £ Cj < 1 



(A.53) 



i>o 



(where x denotes the characteristic function) then we have an absolutely convergent cluster 
expansion given by 

^n^)= ex p{E^y 9 ^K...,w n )n«w}. (a.54) 



Here Wj(S) is the space of circuits compatible with the ion configuration S and having 
j jumps, and (p^i^i, • • • , oj n ) is a combinatoric function on families of circuits whose value is 
zero whenever . . . ,u n } is not a cluster (i.e., whenever the support of the set of circuits 
{ui, . . . , uj n } is not connected in R d+1 ). 

Furthermore, we shall use the following bound, ( A.55|) , which is an extension of Lemma 
3.5 of 



n n n It lb 

9 y ■■■ y x[(x,T)eu 1 ]\( P l(ui,...,u n )\l[\z{u k )\<(n-l)\]lC. 

u>i£W h (S) u> 2 £W j2 (S) u n eW jn (S) 



k=l 



k=l 



(A.55) 



We show that (|A~53D holds in our case as follows. For j ^ 0, (i.e., j > 2) we have 
^ \z{u)\e\ u \ X [{x,r)eu) < ^ t j e~ (c/ ~ 1)|a;| x[(x, r) G lo] 



< (2dt) j [ dTe- {u - 1)T ] J 



< 



{2dt) j ( 



U-l 

For U large enough, U — 1 > cU, for some constant c with < c < 1. Hence 

RHS of (A.56) < Cj 



(A.56) 



(A.57) 



Rigidity of FK Interfaces 



64 



where we define 

Cj := Q) J for j > 2. (A.58) 

For j = we find 

C = e~ PcU (A.59) 

From (|A.58 ) and ( |A.59| ) it follows that the bound (|A.53 ) is satisfied for U and (3 large enough. 
Hence we obtain an expression for the effective hamiltonian: 

1 1 r 11 

H eS (P,S) = -U\A\--^2- j£ <f%(ui,...,u n )Y[zM. (A.60) 



k—~\ 



We define the support, suppa;, of a circuit co, to be its orthogonal projection onto the 
plane r = 0. Hence suppcu C A. 

Let A C A. The potential $a(P, Sa) is introduced as follows: 

-y -y p n n 

MASa) = -oEi 4 x[[J sw P"P UJ k^A](pl(u 1 ,...,u n )l[z(u k ), (A.61) 

P »>1 Wilt*, fc = 1 fe=1 

where 5,4 denotes the restriction of the ion configuration S to A C A. 

Lemma 19 There exists positive constants Uq » t and /3q, such that for all U > Uq and 

(3>(3 

1^x09, 5a)| < Pi ^q=r /H^l>2, (A.62) 

for some constants C\ and c 2 , with n(A) being the minimum length of a closed path which 
passes through all sites of A. 

Proof: Let ji, ■ ■ ■ ,j n denote the number of jumps for the circuits oux, . . . ,u n , such that 
{u>i, . . . , u) n } forms a cluster with support equal to A. 

If Ufc=i supp Uk = A then J2k=i 3 k > n(A) for any u)i, . . . , ui n , then 

IW.JI £ feh £ >o 4 - ^ 

" n>l 31, — ,Jn>0 . . . . 

^ fcJfe >n(A)^ ieW/ w( 5 ) ^eW Jn (S) 

n n 

XtU SU PP^ = ^]IVn( w l>--->^n)| II l*( W k)|- ( A - 63 ) 
k=l k=l 

Now, 

/•/3 1 

22 dr x[(x,r) E uj^- — 1=1, 
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for any uj\ with supp uj\ C A. Introducing this identity into the above equation (|A.63| ) yields 
\*aW,S a )\ < ^EE3 E [dr 4 xKx^eu,]^- 

x ^ ... jf |^(o;i,...,w n )in|z(w fc )|. (A.64) 



Note that |A| > 2 implies that jj > 2 for at least one % 6 {1 . . . n}. We have assumed that 
Ji > 2. 

Let us define z'(u,U) to be equal to the weight z(u), but with explicit dependence on 
the coupling constant U. Let us define 



F(U u ...,U n ) := ^ X [(x,t)eu^) 1 ^^ $ ... ^ 

wjeWj^S) w 2 6VV,- 2 (S) Un€W 3n (S) 

n 

x ... , w n ) | n W{wk, Uk) I (A.65) 



fc=2 



and obtain an upper bound for this function. Eventually, we shall set Ui = for alH = 1 . . . n. 
Note that F(oo, . . . , U n ) = 0. Moreover, 

d 



Hence, 

, 9 



-F(U u ...,U n )\ < jf ^ ... ^ 

n 

x |^(wi, . . . , u n ) I n \z'(uj k , U k ) I (A.66) 



fe=i 



From ( |A.55|) we have that 

r/(^i ; --->^nj| < (n-Ij! 

fc=i 



(j?7 /(^ ■■■ ,«/»)!< (n - 1)! II (A.67) 



where Cj fc (L4) is the same as Cj fc , except for the explicit dependence on Uk- Since 

F(U X , . . . , C/ n ) = - J~ -^F(V, . . . , c^)dy 
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we have that 



F(U u ...,U n ) < J™\-^F(V,...,U n )\dV 



< (n 



n poo 

^■UC jk (u k ) c n (v)dv 

k=2 JUl 



From ( A.58|) we have that 



c n (v)dv < (2dt) ji I {cvy^dv 



(since ji > 2). Hence, 



< 



im, . . . , ty i < (» - i)!^[^] J1 "' n [^r n cm. 



2dUh 1 



ji-i' 



u{ 



2dUh-l ™ r2dtij k " 



fc=2 ^fc" 



fe=2 



Let i be the number of circuits without any jump. From ( |A.64| ) we we have that 

\*a(P,s a )\ < E f\F{U x ...U n )\ 



XgAn>l '"■ Jl,.-,Jn>0 



Further, we have that 



^, ,,.,Co E E ( 2 */ c ) 

„>1 n i=0 l n m=n(A) Ji,-,Jn-i> 2 



E i<2 f 

h 3r>2 



2dt 



c?7 



m— ; 



It is true for all j when m — 1, and by induction 

]_ _ j < 2 m ~ r 2~-' ? '+ :L = 2 m ~( r+1 ) 

j 1 ,... Jr + 1 >2 j r + i >2 Jl,-,Jr>2 JV + 1>2 



Using the bound on the RHS of flOID , we obtain 

1 n | oo 



^[E^+aol e 



m=n(A) 



(2dt/c)' 

jjm—l 



n>l 



m>n(A) 



(Adt/c) 



U 



m—l 
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Now, since Cq 



-(3cU 



[ ([A.59| )1, there exists positive constants Uq » t and /5q, such 



that for all /3 > /3q and U > Uq the sums over n and m on RHS of ( |A.73| ) converge. Thus 
we obtain the bound ( |A.62| ) for some positive constants C\ and c 2 . 

Remark. For A C A we have that Q2.10| ) 



= - 1 

Set t = 1. Then the bound ( |A.62|) can be written as 



\*a(P,S a ) 



C 2 \9(A) 



(A. 74) 



(A.75) 



for some positive constants C\ and C2, with C2/U < 1. This is our desired estimate. 



Appendix B: Proof of Lemma |5] 

From the standard results of cluster expansions |3B], |42|, [^, |T3| it follows that a sufficient 



condition for the convergence of the series given in ( 3.42 ) is given by 



52W(D)eW<l, (B.I) 

Hence, the task of proving Lemma [5] amounts to proving that the condition (|B.1D is satisfied. 
In order to do so, we consider an auxiliary polymer system whose elements (the polymers) 
are denoted by q and defined as follows. 

0:=( 7 ,M 7 ) (B.2) 

where 7 is a contour and M 7 is a decoration of 7. A decoration of 7 is a (possibly empty) set 
of bonds {-Bi, . . . , Bj} C B, such that each bond intersects 7, i.e., Bi n 7 7^ for 1 < i < j. 
Hence a polymer q consists of a contour 7 and a finite set of bonds which intersect its 
support. Let £> 7 denote the set of all bonds which intersect the contour 7. Then M 7 C B 1 . 
For a polymer £> = (7, M 7 ), we define 

I<?|:=l7l+ E (B.3) 

The weight of a polymer is given by 

w(q) := e~ m J] (e"^ - 1), (B.4) 

B£ My 

and satisfies the bound |w(^)| < w (g) where 

Wo (g) .= J] (e^ XB(B) - 1). (B.5) 

BeM 7 
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Each decorated contour D can be considered to be the union of a finite number of intersecting 
polymers, i.e., a connected cluster of polymers. The weight of a decorated contour can then 
be expressed in terms of the weights of its constituent polymers. The decomposition of a 
decorated contour into polymers is however not unique. The condition (|B.1| ) for the model 
of decorated contours can be transcribed into a condition for the auxiliary polymer system, 
by making use of the "tree-graph approximation" used in cluster expansions E2L R6[. We 



sketch the idea below, following | 36| . We first bound the LHS of (|B.1|) in terms of a sum 
over polymers: 

1 n 

£ W(D)eW < £ - £ IW^ 1 - (B.6) 

D30 n>l ei,-,en 7=1 

- ei u...ue„=D 

D 90, connected 



The fact that (|B.6|) is not an equality but only a bound is due to the non-uniqueness of 



decomposition of a decorated contour into polymers. A decorated contour D consisting of n 
polymers, Qi, ■ ■ ■ , Q n , can be represented by a connected, oriented graph, whose vertices are 
the polymers and the lines between pairs of vertices corresponding to intersecting polymers. 
The graph is oriented by introducing an ordering of the vertices. Let Q n be the corresponding 
complete graph, i.e., the graph with n vertices, with a line between each pair of vertices. The 



sum over the polymers on the RHS of (|B.6|) can be bounded by a sum over all tree graphs 



of the corresponding complete graph. Hence, we can write 



1 n 

EW D| <Ei E E---EII^>H (B-7) 



D30 n>l n - TcG n ei30 £n j=l 



where T denotes a tree graph. 

From Lemma 3.5 of IBBH it follows that 



where 



E E-Eil4i> w <(n-ir, (B.8) 

Tcg n gi 30 Qn j=l 



C:=E> (£)e 2|el , (B.9) 
e30 



with Wq(q) being defined through ( |B.5| ). Hence, 

Y,W{D)e\ D \ <E— • ( R1 °) 

D30 n>l 71 

The series on the RHS of ( [B10D converges to C/(l - C) if C < 1, which is < 1 if C < 1/2. 
Hence the proof of ( |B . 1| ) , and hence of Lemma 0, is completed by proving the following 
lemma. 
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Lemma 20 For each C > 0, there exist constants Ao and bo such that for all X < Xq and 
(3X > bo one has the bound 

Y,Me)e 2l * l <C, (B.ll) 

£30 

where w (g) is given by ( \B.3j ). 

Proof: For convenience we consider the sum 



Z pol := $> te)e o|el < C', 

£30 



(B.12) 



for a constant a > 0, and set a = 2 at the end of the proof. Let ko be the smallest positive 
integer for which 

C 2 (3(Xc d e a ) k ° < 1, (B.13) 

where c d = 36. We define 

a := C 2 f3(Xc d e a ) k ° (B.14) 

Let x = C2/3X k and Xq = C2/3X k °. Then for k > k we use the bound e x ' — 1 < xe x and for 
k < k we use the bound e x — 1 < e x . It follows from ( |B.3| ) , ( |B.5| ), and the above bounds, 
that 



w (Q)e al < 



< e -^iA]7l e a|7l x 



n 

B6M 7 
g(B)<k 



e /3C 2 \°W e ag(B) 



X 



J] PC 2 X^e^ X3iB) e a ^ 



BeM-f 

g(B)>k 



(B.15) 



We have that 



\9(B) \g(B) 

< \l\±\ 9{B) <\l\j: f A fc <| 7 |E^ 

B fc>3 B: k>3 

- g(B)=k 



^B.16) 



The symbol J2 x ey denotes the sum over all sites x G A for which at least one nearest 
neighbour bond of the lattice, which contains the site x, is intersected by a face in 7. There 
are 2 1 | such sites in the A. The symbol J2*b denotes the sum over all i?'s containing a fixed 



point. For the last inequality we used the Konigsberg Bridge lemma pp 464-465], which 
gives Cd = (2d) 2 , d being the dimension of the lattice, i.e., c d = 36. 
Hence, we obtain the uniform bound 



n < 

BeM, 



(1C 2 \ 9{B) < e /3C 3 A 3 | 7 | 



(B.17) 
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with 

provided 

For convenience we define 



C s := 



1 - c d \ 
c d \<l. 



a := /3(dA - C 3 A 3 ) 
Then, if ( |B.19| ) is satisfied, we have that 



(B.18) 

(B.19) 
(B.20) 



w { 6 )e a \> 



< 



e -aoh\ x 



n 

seM 7 

s(B)<k 



^9{B) 



1] /3C 2 A s ^e^ 2AS(B) e a ^ B ) 



seAf 7 

9 (S)>fe 



(B.21J 



fB.22) 



Hence, 



emq)^ < e4 e e_aol7l [ n 

£)90 i>0^' -y,B 1 ,...,B j l<i< 3 



j>0J- 7,3!,...,^ 

7USxU...UB :J -30 

[ n ^ 2 A^ fi ) e 

l<i<j 
9(S)>fe 



l<i<j 
ff (S)<fc 



(B.23) 



The sum on the LHS of ( |B.23| ) is over all polymers which contain the origin. For each term 
in the sum on the RHS of (B.23j ), the origin can be contained either in the contour of a 
polymer and/or in one or more of the bonds intersecting it. This sum can be bound by 
a sum over all polymers for which the origin is contained in their respective contours. We 
refer to such a polymer as a pinned polymer. More precisely, a pinned polymer g is defined 
by the connected sequence 7, B x , . . . ,Bj such that 7 3 0. The contribution of each pinned 
polymer must be multiplied by the number of translations containing the origin. Since for 
all i, Bi fl 7 7^ 0, this number is bounded by 



Bi \ 7 | < ]T \Bi\ - 1 < ]T (g{Bi 

l<i<j l<i<j 



(B.24) 



since \B\ < g{B). The factors of "—1" arise from the fact that each bond B must intersect 
the contour 7. Starting from a given contour 7 containing the origin, a pinned polymer 
can be constructed by successively adding bonds which intersect it. For a given 7, the sum 
over decorations can be split into a sum over decorations with 'small' bonds, i.e., B with 
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g(B) < ko, and a sum over decorations with 'large' bonds, i.e., B with g(B) > ko- This gives 
the following estimate: 



-ao|7| 



E 



790 j>0,{S 1 ,...,S 3 }CB 7 j'>0,{Bj S'JCS^ 1<»<J 



E [ E - 1) + E isW - 1; 



i<»<j' 



9(Bi)<*0 



s(s')>fc 



n 



o a a( B i) 



(B.25) 
(B.26) 



n /^a^w^ 

l<i<j J |_l<i<j' 

If both j > 1 and j' > 1 then (since g'(-B) > 3) the sum of the two sums in paranthesis 
in ( p.25|) can be bounded by the product of the two sums, such that the first sum can be 
combined with the first factor of ( |B.26[ ) and the second sum with the second factor. The 
result is 



-oo 7 



1 + 



E (Es(A)-i 



730 



X < 



i>i,{fli Bj}cB 7 l<i<j 

9(Bi)<k 



n 

i<i<j 



o a 9{Bi) 



(B.27) 



E ( E 9W - 1 

g(B',)>k 



j'>l,{B' 1 ,...,B'. / }cB 1 l<i<j' 



i<*<i' 



(B.28) 

The cases with j = or j' = have been incorporated by adding 1 to each factor. Next, we 
estimate the sums in ( |B.27| ). 

The sum in (|B . 2 7| ) can be treated with a 'reverse' high-temperature expansion, i.e., 



resummation, as follows. 



E ( E 9(3,) - i; 

J>l.{Bl B J }CB 7 l<i<j 

9(Bi)<k 



n < 

i<i<j 



ag(Bi) 



< E ( E 9(Bi)) 



j>l,{B 1 ,...,B j }GB 1 l<i<j 
9(Bi)<k 



n 



o a g( B i) 



l<i<j 



d_ 

da 
d 

da 



E 

3*>0,{Bi,...,B i }cB 7 
3(B I )<fc 



n ' 

i<i<j 



ag(Bi) 



H (e a9{B) + l) 



BeB~f 

g(B)<k 



(B.29) 
(B.30) 
(B.31) 

(B.32) 
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pO-g(B) 

E 9 ^e^WTT n < 

g(B)<k g(B')<k 



1) 



< f e n 

BGB 7 B6B 7 
9 (S)<fc 9(S)<fco 

= [ y: g(B)}M a ' e ^) 



(B.33) 
(B.34) 
(B.35) 



ses 7 

9 (S)<fc 



seB 7 

ff(B)<fc 



with a' = a + (log2)/3, guaranteeing that for all g(-B) > 3, exp(ag(B)) + 1 < exp(a'g(B)). 
The estimate of the sum over 'small' decorations can now be completed by using the bound 



n(B) h ° * h ° 

E g(B) = E E nV- |7lE E fc <NE^<^l 



ses 7 

g(B)<k 

with 

At this point we have: 

z pol < E e 

730 



££7 B3i 

9 (S)<fe 



fc=3 s 
g{B)=k 



k=3 



-a \-y\ 



C 4 = (k + l)(c d ) k0+1 . 
l + C 4 \~f\e cW ^' 



(B.36) 

(B.37) 
(B.38) 



x 



1 + 



\ 



E ( E 9{B[) - 1) II /SCalT^Je^ 

j'>i,{sJ,...,b'.,}cb 1 1<»<j' l<i<i' 

V 9(S^)>fc 



(B.39) 



The second term in ( P-39|) represents a sum over the 'large' decorations. Let us denote 
this term by L, i.e., 



L := E ( E ffM - 1) n /W-'^W. 
y>i,{Bj b',}cb 7 l<i<j' i<i<i' 

9(B')>fc 



(B.40) 



and find an upper bound for it. We first make some simplifying estimates. Using ko > 3 and 
I7I > 6 one can see that 

1 + C 4 |7|e C4a ' 171 < C^\e Cia "^ (B.41) 



with, 

a" = a + 1/4. 
For convenience we introduce the notations 



'7>>*:o 



{BeB I Bnj^Q),g(B) > k Q } 



(B.42) 
(B.43) 
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and 

Using the trivial observation 



(3 = (3C 2 , z = Xe a 



E (0<^ E (•) 



we obtain 



{B 1 ,...,B J }cB- h>k() B 1 ,...,B j eB 7 , >kQ 



l < E^ E ( E m) - 1) n ^ m 

j>l J- B 1 ,...,B j eB^ !>ko l<i<j l<i<j 

< E E ••• E tn ^ (fl() ] 
E ( E - i)Pz 9{Bj) - 



X 



We have that 



BjGB 7i>fco l<i<j 



E ( E g(Bi) - i)pz«™ < E E ( E - m 



X£7 SjSa: l<i<j 

s(s,)>fco 



BjeB J:>ko l<i<j 

The sum over .Bj is independent of the choice of the point x. Hence 



RHSof(B.47) < 2\-y\ [g(B 1 )-l + ... + g(B j . 1 )-l + k j -l] 

kj >ko 



(B.44) 



(B.45) 



(B.46) 



(B.47) 



(B.48) 



Iterating the above steps, we obtain the following upper bound to the contribution of the 
large decorations: 



L<J2 l -f E ••• E 

j>l •> ' fei>fco kj>k{) 



J * 



(E^-^IK E 0**) 



i=i 



= 1 Bi 

a(Bi)=ki 



(B.49) 



Since, ki > k > 3, we have that 



j 

I 

i=i 



E(h - 1) < U( k * - V 



(B.50) 



j-i 
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Each factor of (k{ — 1) can be inserted into the sum over fcj. This yields 



\3 3 



L < E^-IKE(^-i) E P* ki ) 

j>l J" *=1 fci>fc 



g(S i ) = fc i 



j>l J- i=lki>k 



< E^f^o(c^) fco 
i>i J ' 

< E^ofe^] 



E 

fe>i 

(1 - c d zf 



provided 

c d z = c d e a \ < 1 
By definition (( pl3D , (gig) ) we have that 

/3(c d z) fco = /3C 2 (c d e a A) fef > = « < 1 

Hence, i/ QB.52|) holds, then 



< exp 



\l\k c d z 



L - E ;. \ , , ,. ,2 



(1 " c d zf 
exp(ai|7|) 



with 



ai := 



k c d z 



(1 - Qz) 2 

Hence from ( [Egg ), ([BHD and (|B34| it follows that 



^poi < 5>- ao|7l C 4 |7|e C4a " 171 [l + e ai171 

790 

< ^ e - aol7l 2(l74|7|e C4a " l7l e ai171 , 



790 



since a\ > 0. Defining 



g := a - log c d - ai - a"C 4 



(B.51) 
(B.52) 
(B.53) 



(B.54) 
(B.55) 



(B.56) 



(B.57) 
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we have that 



fc>6 

e~ c - 

< 2C 47 



kq 



-q\2 



(B.58) 
(B.59) 



provided 

q>0. (B.60) 

Let us inspect the condition ( P.60|) in more detail. From the definitions, (B.20| ), ( |B.55|) , 
flggg) and flET37D , of a , ai, a" and C 4 , and (|B~oD|) , it follows that (|B39]) holds provided 



A/3 X(A) > a + log c d + k 



c d Xe a 



where 



X{\) := C 



c d Xe a ) 2 
C2C d \ 2 



+ (a + -) (k + l)cg° +1 , 



(B.61) 



^B.62) 



1 - qA 

and k satisfies the bound ( [B. 13| ) . Before proceeding further, let us recall the conditions 
which have been imposed on A in order to arrive at the form (|B.59|) . These are given 
by ( |B~T9D and flB32l ). Moreover, for flBUTD to hold, it is necessary that X(X) > 0. The 
inequalities (|B.19| ) and (p.52| ) are satisfied if 



A < Ai :-- 



[c d e 



a\-l 



[B.63) 



Since X(X) is monotone decreasing, X(\) > can be satisfied by requiring A to be smaller 
than the solution of X(X) = 0, i.e., 



A < A 2 := 



c d 1 



In order to satisfy ( B.63|) and ( p.64|) we choose A > Ao with 

Ao^OB^eT 1 

where 



B :-- 



1 + 



4c d C 2 



C-i 



(B.64) 

(B.65) 
(B.66) 



Note that B > 1 since C\, C2 and c d are positive. 

The quantity q defined through ( |B.57[ ) can be expressed as a function of A and b := (3X. 
Next, let us determine bo such that for all b > bo and A < Ao, ( B.60|) is satisfied. On the 
RHS of ( |B.61p , k Q (expressed as a function of A and b) has to satisfy 



k = k (X,b) > 1 



log (C 2 c d e a b) 
\og(Xc d e a ) 



[B.67) 
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This follows from the defining relation ( B.13 ) for ko- For all A < Ao this can be achieved by 
choosing 



ko = k {b) 



1 - 
1 + 



log (C 2 c d e a b) 
log (A c d e a ) 
log (C 2 c d e a b) 
log(5c d ) 



^B.68) 



where B is defined through ( |B.66| ). Let us denote the RHS of ( |B.61[ ), with the above choice 
k , by A(\,b), i.e., 



k Xe a 1 — — 

A{X, b) = a + log c d + _\ dXe ay + c ^ a + 4)^° + l )° k d- 



(B.69) 



Recall that X(X) is monotone decreasing, and note that A(\,b), for a fixed value of b, is 
monotone increasing in A. Hence, for all A < Aq we can satisfy ( |B.61|) by requiring that 



bX(X )>A(X ,b). 



(B.70) 



That ( |B.70|) is satisfied for all b > bo, for some 60 > 0, follows from the fact that A(X Q ,b) 
increases strictly less than linearly in b. More precisely, 



A(X , b) = A + Ailog (d b) + A 2 (d b) r + A 3 (d b) r log (do b), 
where Aq, Ax, A 2 , A 3 and do depend on C\, C 2 , e a and c d , and r is given by 

logc d 



In particular, 



Since B > 1 (see 



log Bc d 



do = C 2 c d e a 



(B.71) 



(B.72) 



(B.73) 



we have that r < 1. This proves the convergence of the series for 
Z pol and the bound ( p.59|) . 

It is now easy to see that, in fact, the bound C for Z pol as in the statement of the lemma, 
can be made arbitrarily small by choosing 6 arbitrarily large. With k = k Q , C is given by 



C'(b) = 2c 2 d 2 + 



\og(Bc d ) 



(Cxb) r x 



-q(b) 



(1 _ e -<?(fe))2- 



(B.74) 



where we have used the definitions of C 4 [( [B.37| )1 and k [( |B.68| )1. The constants d and B 
are defined in ( |B.73|) and (|B.66|) . The function q(b) is bounded below by bX(X ) — A(X , b), 
which increases linearly in b. Then again since r < 1, C'(b) — ► as b — > 00. | 
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